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Introduction 



Quantum spin chains are probably the simplest quantum mechanical systems 
showing a wide variety of interesting properties [H [2l [3l |4l [5] , a main one being 
the existence of quantum phase transitions, that is transitions at zero temper- 
ature driven by large quantum fluctuations [B]. Singularities occuring strictly 
[7] at the zero-temperature transition point can, however, produce a typical 
signature at (very) low temperatures. On the experimental side, physicists are 
nowadays able to produce artificial samples with behaviour that fits very well 
with the theoretical descriptions [H HI [10] . 

The low-dimensionality of these systems allows the use of efficient analytical 
and numerical tools such as the Bethe ansatz [TTl [12] , bosonization [13l [14] , 
or fermionization [151 [T] , exact diagonalization of finite chains or the numerical 
density matrix (DMRG) approach [T71 [TS] . We focus our attention in this 
review on free-fermionic quantum spin chains, which are spin models that can 
be mapped on systems of non-interacting fermions [T]. The free nature of the 
elementary excitations allows for an exact diagonalization of the Hamiltonian. 
It is not here necessary to argue on the usefulness of exact solutions for the im- 
plementation of the general comprehension we have of such many-body systems. 
In particular, they can be used as a garde fou when dealing with more complex 
systems, unsolved or unsolvable. Moreover, some of the properties they show 
can still be present in truly interacting systems, that is, models that cannot be 
or have not yet been mapped on free particles or solvable problems. 

Most of the studies during the last decades, have been devoted to understand 
the influence of inhomogeneities [19] , such as aperiodic modulation of the cou- 
plings between spins |20j or the presence of quenched disorder [21 , on the nature 
of the phase transition [52]. This has culminated in the work of D. Fisher [23] 
on the random transverse field Ising quantum chain. The extremely broad dis- 
tribution of energy scales near the critical point, for such random chains, allows 
the use of a decimation-like renormalization group transformation [24j . Another 
aspect which is considered in this paper is the non-equilibrium behaviour of such 
fermionic spin chains. The focus will be on homogeneous systems, since they 
allow analytical calculations, although expressions are given which are valid for 
general coupling distributions. 

The paper is organised as follows: in the two next chapters, we present the 
general features of the free fermionic spin models and the canonical diagonal- 
ization procedure first introduced by Lieb et al. [1 . A detailed discussion is 
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given on the excitation spectrum and the associated eigenvectors. It is also 
pointed out how one can extract the full phase diagram of such spin chains 
from the knowledge of the surface magnetization. Results on dynamical cor- 
relation functions are reviewed. Chapter [3] gives some rapid introduction on 
the studies done on the aperiodic Ising quantum chain. We have focused the 
attention on the anisotropic scaling and the weak universality found in such 
systems. In chapter 21 we present results obtained by ourself and others within 
the random Ising quantum chain problem. After quickly reviewing RG-results 
worked out on normal homogeneous distributions of disorder, we present the 
results obtained on the critical behaviour of random chains with either Levy 
type disorder or inhomogeneous disorder. The attention is concentrated on the 
surface critical behaviour since, due to the particularly simple expression of the 
surface magnetization, it is possible to obtain many exact results for that quan- 
tity. The following chapter deals with the non-equilibrium behaviour. After 
solving the Heisenberg equations of motion for the basic dynamical variables, 
we present some aspects of the relaxation of the transverse magnetization. We 
show that systems with either conserved or non-conserved dynamics present, 
however, some similarities in their relaxation behaviour. This is illustrated 
with the XX-chain for conserved dynamics, and with the Ising chain for non- 
conserved dynamics. Two-time functions are also considered and aging, that 
is, the dependence of the relaxation process on the age of the system, is also 
discussed. 



Chapter 1 

Ising Quantum chain 



1.1 Free fermionic models 



The generic XY Hamiltonian that we will consider is UHl [57] 



H 



L-l 



1 + K 
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(1.1) 



n=l ^ 



where the 



= l(g)l---(g)cr(g)l---(g)l 



are Pauli matrices at site n and n is an anisotropy parameter with limiting 
values K = \ corresponding to the Ising case with a Z2 symmetry and k — Q 
describing the XX -model which has U{1) symmetry. We will consider here only 
free boundary conditions. The phase diagram and the critical behaviour of this 
model are known exactly since the work of Barouch and McCoy in 1971 |28l[29| 
who generalized results obtained previously in the case of a vanishing transverse 
field T] , or at K = 1 [27] . It was first considered in the framework of conformal 
invariance [30j in Ref. [31j . In this section, to give a self-consistent presentation, 
we present in full details the diagonalization procedure, following, more or less 
closely, the initial work of Lieb et al. pj . 

The Hamiltonian can be mapped exactly on a free fermion model, consisting 
of an assembly of non interacting Fermi-Dirac oscillators. To proceed, let us 
first introduce the ladder operators 



In the diagonal basis of the cr^ component, they simply act as cr"*"! |) = | t) and 
a~\ T) = I i)- They satisfy the anticommutation rules at same site 



a 



-(a-±ia«). 



a-} = l 
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and by construction they commute at different sites. They look like Fermi 
operators apart that they commute on different sites. True fermionic operators 
are obtained through a Jordan- Wigner transformation [15] . 

Consider a basis vector with all spins down (in z-direction) and use the 
notations |1) = | t) and of course |0) = | i), then this state will be the vacuum 
state destroyed by all the lowering operators cr^ : 

cr-|00...0> ==0 Vn. 

From this vacuum state, all the other states can be built up by applying raising 
operators cr+. The situation looks very much the same than with fermions. But 
to have really fermions we need antisymmetry, i. e. anticommutation rules also 
at different sites. We introduce new operators 

c„ = A{n)a~ , 
c+^a+A+{n) 

where A{n) is a unitary operator commuting with : 

[A(n),a±]=0. 

Then automatically we have 

{c+,c„} = {o-+,cr,^} = 1 

and in particular 

< = 2a+<j- - 1 = 2c+c„ - 1 . (1.2) 

In order to fulfil the antisymmetry principle (sign change under particle ex- 
change), the creation and annihilation operators must satisfy ^32j 

c^|ni,n2, ...,n;, ...) = (-l)^'(l - n()|ni, ^2, n; -|- 1,...) , 
c;|ni,7i2, ...,n;, ...)(-l)^'n;|ni,?i2, ...,ni - 1, ...) , 

where — is the number of particles on the left of the / site. We should 

find now an operator representation of the sign factor (—1)^'. It is obvious that 
the choice 

m = n(-<) 

i=l 

will perfectly do the job I32i- This leads finally to the so called Jordan- Wigner 
transformations JS] 

1 n—1 

1 = 1 4=1 

and the adjoint relation. Inverting these relations, we obtain expressions of the 
ladder operators and the original Pauli matrices in terms of the fermionic cre- 
ation and annihilation operators c+ and c. Replacing this into our Hamiltonian 
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we obtain the quadratic form 

n,m— 1 

where the real matrices A and B are respectively symmetric and antisymmetric 
since the Hamiltonian is hcrmitian. The quadratic nature of the Hamiltonian 
in terms of the Fermi operators insures the integrability of the model. As in 
the bosonic case, this form can be diagonalizcd by a Bogoliubov transformation 

mm- 

Let us mention here briefly that more general quantum spin-1/2 chains are 
not tractable using this approach. For example, in the Heisenberg model [11] 

H ~ J ^ fT„(T„+l 
n 

there are terms of the form 

z z + + 

leading to interacting fermions. With longer range interactions, such as a^a^^p, 
or magnetic fields in the x or y directions, another problem arises due to the 
non- locality of the Jordan- Wigner transformations. For example a next nearest 
neighbour interaction cr^cr^_|_2 generates quartic terms like 

1 Ct! + 1 C„_|_ 2 

where the fi and v upper-scripts refer to either creation or annihilation operators. 
In the case of a magnetic field, let us say in the x direction, we have additional 
terms proportional to the spin operators 

< = + ^^)(^^ - ^^)) + 

which are clearly even worse. Nonlocal effects also appear when dealing with 
closed boundary conditions [U [26] . 

Nevertheless, the Hamiltonian (|3.5p is not the most general free-fermionic 
expression. We can still add for example terms of the form a^a^^^—a^a^^j^, the 
so called Dzyaloshinskii-Moriya interaction [SH [35] or play with the boundary 
conditions [36l. 



1.2 Canonical diagonalization 

We now come to the diagonalization of Hamiltonian (|3.5p . For that purpose, 
we express the Jordan- Wigner transformation in terms of Clifford operators [37j 
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= n 

= (1-5) 

These operators are the 2L-generators of a Clifford algebra since 

{r;,ri.} = 25y J„fe , Vz,j = l,2;Vn,fc = l,...,L (1.6) 

and are real operators, that is r*j^ = r*j. They can be viewed as non-properly 
normalised Major ana fermions. The different terms in the original Hamiltonian 
are expressed as 



ir^r^ 



l-L n+1 



Introducing the two-component spinor 

r„ = 



^ 71 



(1.7) 



one can write the Hamiltonian in the form 

L-l 

4 



n=l 



where = (r,\,r2 ) and 



(1.8) 



-i\ 




^ 




i i 






I) 



are the Pauli matrices, not to be confused with the initial spin operators. In- 
troducing the 2L-component operator = {t\,t\...,v\)^ the Hamiltonian is 
given by 

H = -T^TT 

4 



where T is a 2L x 2L hermitian matrix 



/ D F 

Ft D 

'■■ 





V 




F 



given by 

\ 
... 



Ft D F 

Ft D y 



(1.9) 



(1.10) 
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with 



n = hay , F = - (cry + ikct^) 



(1.11) 



To diagonalize H, we introduce the unitary transformation matrix U build 
up on the eigenvectors of the T matrix: 

TV,^e,Vg ,q^l,...,2L 

with the orthogonahty and completeness relations 

2L 2L 

J2 V* mg, (*) = V , E K (*)^?(*') = '5-' • 

1=1 q=l 

Inserting into (|3.12p the expression T = UAU^^ where \pq — e^Spq is the diag- 
onal matrix, one arrives at 

1 1 \ 

H = -rtUAUtr = -Xt AX = - ^ e^x+x^ 

9=1 

with the diagonal 2L-component operator 

X = u^^r . 

We introduce now the following parametrisation, which will become clear later, 
of the eigenvectors Vq-. 

( U^) \ 

0,(2) 
-iV.,(2) 



V2 



(1.12) 



Utilising this parametrisation, the operators Xq and their adjoints are given by 



1 ^ 



Now, if we consider the normalised operators 
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we can easily check that together with the adjoints they define Dirac 
Fermions, that is they satisfy the anticommutation rules 



{d+,dg,}^6,^,, , {d+,d+}=0, {dg,dg,} = 0. 

For example, the first bracket is evaluated as 



(1.13) 



- i^,(*)0;,(j){r2,r]} + i0,(z)v*,(j){r.^r2} 

and using the anticommutation rules for the Clifford operators and the normal- 
isation of the eigenvectors one is led to the above mentioned result. 
Finally we have the free fermion Hamiltonian 



^ 2L 



(1.14) 



9=1 



We now take into account the particular structure of the T matrix. Due 
to the absence of the diagonal Pauli matrix cr^ in the expression of T, the 
non- vanishing elements are those with i + j odd, all even terms are vanish- 
ing. This means that by squaring the T matrix we can decouple the original 
2L-eigenproblem into two L-eigenproblems. The easiest way to see this is to 
rearrange the matrix T in the form 



C 

ct 



(1.15) 



where the L x L matrix C is given by 



( h Jy 

Jx h 

J X 

O 



J. O 



J II 



(1.16) 



V Jx h J 

with Jx = (1 + k)/2 and Jy ~ (1 — k)/2. Squaring the T matrix give^ 



CCt 

etc 



(1.17) 



^The supcrsymmetric structure appearing in the matrix has been used in Ref. |40| . 
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In this new basis the eigenvectors Vq are simply given by 

0,(1) \ 
0,(2) 

ML) 

V -^ML) J 

and together with we finally obtain the decoupled eigenvalue equations 

CCt0, = e^^q (1.18) 

and 

CtC^, = e^t^g . (1.19) 

Since the CC^ and C^C are real symmetric matrices, their eigenvectors can 
be chosen real and they satisfy completeness and orthogonality relations. This 
justifies the initial parametrisation of the vectors Vq and one recovers the original 
formulation of Lieb, Schultz and Mattis [1]. 

Finally, one can notice another interesting property of the T matrix which 
is related to the particle-hole symmetry [38]. Due to the off-diagonal structure 
of T, we have 

- iClpq = eq(l)q , 
CUq = -ieq^q (1.20) 

and we see that these equations are invariant under the simultaneous change 
eq — !■ —eq and ijjq — > —^^q- So, to each eigenvalue > associated to the vector 
Vq corresponds an eigenvalue eq> = —£q associated to the vector 

M^) \ 

cj,q{L) . 

i^9(l) 

V '^i>qiL) ] 

Let us classify the eigenvalues such as 

Eg+L = -e, Vq = l,...,L 
with > Vg = 1, L. Then the Hamiltonian can be written as 

1 

^ = 2 H (^-J^^? " iqd'^+Ldq+h) 

q=l 



1 
3 
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where the operators with q = 1, ...,L are associated to particles and the oper- 
ators with q = L + 1,...,2L are associated with holes, that is negative energy 
particles. So that, by the usual substitution 



vt=dg Vq = l,...,L 
+ ^dg+L Vq==l,...,L 



we rewrite now the Hamiltonian in the form 



1 ^ 



9=1 



9=1 



(1.21) 



(1.22) 



1.3 Excitation spectrum and eigenvectors 



The problem now essentially resides in solving the two linear coupled equations 
— iCV'ij = iq4>q and iC^^q — eqipq. We will present here the solutions of two 
particular cases, namely the XF-chain without field [Tl[26] and the Ising quan- 
tum chain in a transverse field [27] . In the following we assume, without loss of 
generality, that the system size L is even number and k > 0. 

1.3.1 Xr-chain 

From iC^0q — Eqipq, we have the bulk equations 



i_^0^(2fc-l) + i±^0,(2fc+l) 
^~''-0,(2fc) + i±^0,(2fc + 2) 



-eq%l)q{2k) 

-eqijq{2k + l) (1.23) 



2 ' 2 

Due to the parity coupling of these equations, we have two types of solutions: 
^^(2fc) = Vg(2fc-l) = yk (1.24) 

and 

(l>l\2k - I) = ij'q\2k) = yk (1.25) 
In the first case, the bulk equations that remain to be solved are 
1 - K 



2 yqi2k-l) + ^cPl{2k + l) = -eqi:!^{2k) 
with the boundary conditions 

0^(L+1)=V4(O) = O. (1.26) 
Here we absorb the minus sign in equation (|1.23p into the redefinition 
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Using the ansatz 0^(2fc- 1) = e'^^^fc-i) and ipl{2k) = e^'^^'^e*^" to solve the bulk 
equations (|1.23[) . we obtain 

cosq + iKsing = e^e*^" (1-27) 

that is 

€q = \J cos2 g + k2 sin^ q > (1.28) 

and the phase shift 

9q — arctan(Ktanq) , (1-29) 

with < 6*^ < 27r to avoid ambiguity. The eigenvectors associated to the positive 
excitations satisfying the boundary equations are then 

<^j{2k + l) = A,sin{qi2k + l)-9,) 

^l{2k) ^ A,sin{q2k) , ^ > 

with 

qiL + l)=mr + 9g (1.31) 

or more explicitly 

q= [ n H — arctan (^tanq) I . (1-32) 

L + I \ n J 

The normalisation constant Aq is easy to evaluate and is actually dependent on 
q. Using 9q — q{L + 1) + titt, one can write 0^ in the form 

(t)[{2k + 1) = -~Aq5q smq{L - 2k) , 

where 5q = (—1)" is given by the sign of cosq{L + 1). The equation (|1.32p has 
L/2 — 1 real solutions, that in the lowest order in 1/L are given by 

7r 

Qn^ J^{n- Vn) (1.33) 

with 

n 1 ( (n-KW L 

Vn = T arctan I K tan I — , n = 1, 2, — — 1 . (1-34) 

There is also a complex root of (|1.32p 

TT 

q^^-^w (1.35) 

where v is the solution of 

tanhw = Ktanh[i;(_L + 1)] . (1.36) 
With the parametrisation x = e^^" and = ji^, one is led to the equation 

^2 



1 — a;-^(l — p^x) 
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and the first nontrivial approximation leads to 

The excitation associated with this localised mode (see the form of (pq^ and tpq^ 
with go = 7r/2 + iv) is exponentially close to the ground state, that is 

~ (1 + p2)^L ^ (138) 

From this observation, together with some weak assumptions, the complete 
phase diagram of the system can be obtained. We will discuss this point later. 

The solutions of the second type satisfy the same bulk equations but the 
difference lies in the boundary conditions <pg^{0) — i^g^ {L + 1) = 0. The eigen- 
vectors are given by 



0"(2fc) = Agsm{q2k) 



with 



?/'"(2fc + 1)^ Aq sin(g(2fc + 1) + Oq) = -AgSq smq{L - 2k) 



(1.39) 



q = ^ ^ ( " arctan(Ktan(7) ) (1-40) 



which has L/2 real roots. To the leading order, one gets 



qn = Y^{n- Vn) (1-41) 



n 1 / / riTT \\ ^ ^ L 



with 

Vn = \ — arctan(Ktan( — )) ,n = l,2,...,— (1-42) 

which completes the solution of the XF-chain. 

1.3.2 Ising-chain 

The solution of the Ising chain {k — 1) proceeds along the same lines [2J]. The 
bulk equations are 

h(j)g{k) + (jjqik + 1) ^ eq4>g{k) (1.43) 

with the boundary conditions 

cPq{L+l)^,Pq{0)=0. (1.44) 

With the same ansatz as before, one arrives at 

/i + = e^e''^' (1.45) 

that is 



^J{h + cosq)2 + sin^ q (1-46) 
and 

/ sing \ / 
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Taking into account the boundary conditions, the solutions are readily expressed 
as 

(pq{k) = Asin{qk — 6q) 
ipq{k) — ~Asm{qk) 

where g is a solution of the equation 



(1.48) 



TT / 1 / sino , , / .„s 

q = n + - arctan — . (1.49) 

L + 1 \ TT \h + cosqJ ' 



The eigenvectors can then also be written in the form 

(l)q{k) = -A(-l)" sm{q{L + l-k)) 
ipq{k) = —Asm{qk) 



(1.50) 



In the thermodynamic limit, L — > oo, for h > 1, the equation (jl.49p gives rise 
to L real roots. On the other hand, for ft, < 1, there is also one complex root 
qo — TT + iv associated to a localised mode such that v is solution of 

tanh(i;(£ + 1)) = - , ""^^'^ . (1.51) 
h — cosh V 

To the leading order, we have v ~ \nh. The eigenvectors associated to this 
localised mode are 



(jjg^ (k) = A{-1)'' smh{v{L +l-k)) 
^q^{k) = -A{-lfsin\i{vk) . 



(1.52) 



Exactly at the critical value h — 1^ we have 6q = q/2, which gives a simple 
quantisation condition: 

9= 1,2,...,L (1.53) 

Changing q into tt — g, we have [31] 



{k)^^,{-lf+'cos{q[k-l/2)) 



2 

V2L + 



i^,ik) = ^^,i-irsm{qk) 



and 

9 



eg = 2 



sm — 
2 



(1.54) 
(1.55) 



with q={2n+ 1)tt/{2L + 1) and n = 0, 1, L-l. 
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Chapter 2 

Equilibrium behaviour 



2.1 Critical behaviour 

From the knowledge of the eigenvectors (j) and ip, and the corresponding one- 
particle excitations, we can in principle calculate all the physical quantities, such 
as magnetization, energy density or correlation functions. However, they are, 
in general, complicated many-particles expectation values due to the non-local 
expression of the spin operators in terms of fermions. We will come later to this 
aspect when considering the dynamics. Nevertheless, quantities that can be 
expressed locally in terms of fermion are simple, such as correlations involving 
only operators or af. 

2.1.1 Surface magnetization 

A very simple expression is obtained for the surface magnetization, that is the 
magnetization on the x (or y) direction of the first site. The behaviour of the first 
spin [321 US] gives general informations on the phase diagram of the chain [551 
[29l [30] . Since the expectation value of the magnetization operator in the ground 
state vanishes, we have to find a bias. The usual way will be to apply a magnetic 
field in the desired direction, in order to break the ground state symmetry. Of 
course this procedure has the disadvantage to break the quadratic structure 
of the Hamiltonian. Another route is to extract the magnetization behaviour 
from that of the correlation function. In this respect, the x(y) component of 
the surface magnetization is obtained from the autocorrelation function G(t) = 
(erf (O)cr^(r)) in imaginary time r where erf (r) = e'^^afe~'^^ . Introducing the 
diagonal basis of iJ, we have 

G(t) = |(a|CTf|0)pe-"(^'~^«) +^|(i|af|0)|2e-"(^'-^«) (2.1) 

i>l 

where |0) is the ground state with energy Eo and |cr) = vtl^) the first excited 
state with one diagonal fermion whose energy is Eg- = Eq + ei. From the 



19 



20 



CHAPTER 2. EQUILIBRIUM BEHAVIOUR 



previous section we see that we have a vanishing excitation for /i < 1 in the 
thermodynamic hmit L oo leading to a degenerate ground state. This imphes 
that in the hmit of large r, only the first term in the previous expression of the 
autocorrelation function contributes: 

hm G(r) = K]' , (2.2) 

r — >oo 

where = {a\af\0). Noticing that af — Tl — cf + ci and making use of the 
inverse expression, = J^q 4'q{n){T]q + f]q), one obtains 

m: = (aK|O)=0i(l). (2.3) 

Similarly, one can obtain to| = (CTlcrj^lO) = 

Following Peschel |32], it is now possible to obtain a closed formula for the 
surface magnetization [43j . In the semi- infinite limit L — s- cx), for h < I, the 
first gap — Eq = ei vanishes due to spontaneous symmetry breaking. In this 
case, equations (|1.20|) simplify into 

Ct0i =0 

cv^i = • 

Noticing that changing k into —k, rrv^ and are exchanged, in the following 
we will consider only the x-component. 

To find the eigenvector 0i, we rewrite C^0i =0 in the iterative form 

01 (n) ) " \^ 0i(n- 1) J ^ ' 

where K„ is a 2 x 2 matrix associated to the site n whose expression for homo- 
geneous coupling constants is 

2h 1-n 



K„ = K = -( i+« i+« 1 . (2.6) 



By iterations, we obtain for the (n + l)th component of the eigenvector (pi the 
expression 

0i(n + l) = (-l)>i(l)(K"),, (2.7) 

where the indices 11 stand for the 1,1 component of the matrix K". The nor- 
malisation of the eigenvector, 4>i{i) = li leads to the final expression 



n=l 



l(K")nl 



'1/2 



(2.8) 



First we note that the transition from a paramagnetic to an ordered phase is 
characterised by the divergence of the sum entering (|2.8p . On the other hand, 
since for a one-dimensional quantum system with short-range interactions, the 
surface cannot order by itself, the surface transition is the signal of a transition 
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in the bulk. It means that one can obtain some knowledge of the bulk by 
studying surface quantities. To do so, we first diagonalize the K matrix. The 
eigenvalues are 

^ ± ^/h^ + k2 - 1 



A± = 

1 + K 

for + > 1 and complex conjugates otherwise 

1 



(2.9) 



A. 



1 + K 



h±iVl-h^- k2 pexp(±i79) (2.10) 



with p = ^y {1 — k)/{1 + k) and = arctan(-\/l — h? — n'^/h) and they become 
degenerate on the line h? + n"^ = 1. The leading eigenvalue gives the behaviour 
of 4>i{n) ~ which, for an ordered phase, implies that |A+| < 1. The 

first mode is then localised near the surface. From this condition on A+, we 
see that it corresponds to h < he = I. For h > 1, the surface magnetization 
exactly vanishes and we infer that the bulk is not ordered too. So that for any 
anisotropy n, the critical line is at /i = 1 and in fact it belongs to the 2d-lsmg 
[44] universality class0 As it is shown hereafter, the special line + = 1, 
where the two eigenvalues collapse, separates an ordinary ferromagnetic phase 
(/i^ -(- ^2 ^ -[^-j fj-Qjj^ g^j^ oscillatory one {h'^ + k2 < 1) [28]. This line is known 
as the disorder line of the model. The line at vanishing anisotropy, k = 0, 
is a continuous transition line (with diverging correlation length) called the 
anisotropic transition where the magnetization changes from x to y direction. 

In the ordered phase, the decay of the eigenvector 0i gives the correlation 
length of the system, which is related to the leading eigenvalue by 

0?(n)~|A+p"~exp(^-|) (2.11) 

so that 

By analysing A+ , it is straightforward to see that the correlation length exponent 
defined as ^ ~ ^ with S oc 1 — h for the Ising transition and ^ oc k for the 
anisotropic one, is = 1. Of course, a specific analysis of formula ()2.8|) leads to 
the behaviour of the surface magnetization. At fixed anisotropy k, close enough 
to the Ising transition line we have 

m^,^{l-hf^' (2.13) 

giving the surface critical exponent /3f = 1/2. In the oscillatory phase, a new 
length scale appears given hy 'd~^. A straightforward calculation gives for the 



^ To get an account on the connection between d-dimensional quantum systems and d+ 1- 
dimensional classical systems, on can refer to Kogut's celebrated review 1451 . The idea lies in 
the fact that, in imaginary (Euclidean) time t, the evolution operator e~'^^ , where H is the 
Hamiltonian of the d-dimensional quantum system, can be interpreted as the transfer matrix 
of a d + 1 classical system. 
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eigenvector the expression 



n 



1 + K 

ft, tan -i? 



sin(ni?)0i(l) 



(2.14) 



which leads to 



K 



1/2 



(2.15) 



close to the anisotropy hne, so (3^—1/2 too. 

We have seen here how from the study of surface properties [43] one can 
determine very simply (by the diagonalization of a 2 x 2 matrix) the bulk phase 
diagram and also the exact correlation length [5S]. One may also notice that 
expression (12. 8p is suitable for finite size analysis [15], cut off at some size L. 
In fact for the Ising chain, to be precise, one can work with symmetry breaking 
boundary conditions 47J. That is, working on a finite chain, we fix the spin 
at one end (which is equivalent to set = 0) and evaluate the magnetization 
on the other end. In this case, the surface magnetization is exactly given by 
(|2.8p where the sum is truncated at L — 1. The fixed spin at the end of the 
chain, let us say cr£ = +, leads to an extra Zeeman term since we have now in 
the Hamiltonian the term — "^^^^ o"l-1' where the last coupling Jl~i plays the 
role of a magnetic field. The vanishing oi Hl induces a two-fold degeneracy of 
the Hamiltonian which is due to the exact vanishing of one excitation, say ei. 
This degeneracy simply reflects the fact that [af,H] = 0. On a mathematical 
ground, we can see this from the form of the T matrix which has a bloc-diagonal 
structure with a vanishing 2x2 last bloc. From the vacuum state associated to 
the diagonal fermions, |0), and its degenerate state ?7j^|0), we can form the two 
ground states 



associated respectively to = + and = — . Since we have a boundary 
symmetry breaking field, we can directly calculate the surface magnetization 
from the expectation value of af in the associated ground state. It gives 
(+l'''il+) — 4'i{^)j where (f>i{l) is exactly obtained for any finite size from 
C^4>i — 0. This finite-size expression has been used extensively to study the 
surface properties of several inhomogeneous Ising chains with quenched disor- 
der [ISl Sni HOI [SB , where a suitable mapping to a surviving random walk [JT] 
problem permits to obtain exact results. 

2.1.2 Bulk magnetization 

As stated at the beginning of this section, quantities involving ct^ or cr'' operators 
are much more involved. Nevertheless, thanks to Wick's theorem, they are 
computable in terms of Pfaffians [5^ [55] or determinants whose size is linearly 
increasing with the site index. For example, if one wants to compute (crltrf |0), 
the magnetization at site I, one has to evaluate the expectation value [53' 




(2.16) 



mi = {a\af |0) {Olrj.A^B.A^B^ . . . A, |0) , 



(2.17) 
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where we have defined Ai = Vj and Bi = —iTf in order to absorb unnecessary 
i factors. These notations were initialy introduced by Lieb et aZ [l]. Note that 
B^ — —1. Since Ai and Bi are linear combinations of Fermi operators: 



Vq) 



Sz-^V^W , (2.18) 



we can apply Wick's theorem for fermions J^. The theorem states that we may 
expand the canonical (equilibrium) expectation value, with respect to a bilin- 
ear fermionic Hamiltonian, of a product of operators obeying anticommutation 
rules, in terms of contraction pairs. For example, if we have to evaluate the 
product (C1C2C3C4), we can expand it as 

(C1C2C3C4) = (CiC2)(C3C4) - (CiC3)(C2C4) + (CiC4)(C2C3) . 

Due to the fermionic nature of the operators involved, a minus sign appears 
at each permutation. In our case, it is easy to see that the basic contrac- 
tions (0|v4iv4j|0) and {0\BiBj\Q) are vanishing for i ^ j. The only contributing 
terms are those products involving only pairs of the type (0|77iA|0), (0|77ii3|0) 
or (0|i3A|0). One may also remark that it is unnecessary to evaluate terms 
of the form (0|7/ii?|0) since in this case there is automatically in the product 
a vanishing term (0|^iylj|0) with i 7^ j. The simplest non-vanishing product 
appearing in the Wick expansion is 



{0\ViAi\0){0\BiA2\0)...{0\Bi^iAi\0) . 
The local magnetization m; is then given by the I x I determinant [53l |47] 



mi 



with 



and 



Hi Gil G12 
H2 G21 G22 

Hi Gil G12 



Gi/_i 
G2/-1 



G, 



ii-i 



H, = {Q\mA,\0) = Mj) 



G,k^{0\BkA,\0) = -J2Mj)Mk) ■ 



(2.19) 



(2.20) 



(2.21) 



This expression for the local magnetization enables to compute, at least nu- 
merically, magnetization profiles [Ml [SS] and to extract scaling behaviour. For 
example far in the bulk of the Ising chain we have near the transition the power 
law behaviour mt (1 — ft.)^/^ for /i < 1 and zero otherwise [27J. To conclude 
this section, one can also consider more complicated quantities such as two-point 
correlation functions [Tl I27j in the same spirit. 
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2.2 Time-dependent correlation functions 

Time-dependent correlations functions are of primary importance since exper- 
imentally accessible dynamical quantities are, more or less simply, related to 
them. For general spin chains, the exact analysis of the long time behaviour of 
spin-spin correlations is a particularly difficult task. The generic time-dependent 
spin-spin correlation hmction is 



where i,j are space indices, — x,y,z and where the average ( . ) = 

Tr{ . e~^^} /Tr{e~^^} is the canonical quantum expectation at temperature 
T — 1/(3. The time-dependent operator a^{t) — e'^Vfe"'^' is given by the 
usual Heisenberg representation. Most of the approximation schemes developed 
so far [Sni [57] are not really relevant for such many-body systems, at least at 
finite temperature. One is ultimately forced to go on numerical analyses, basi- 
cally by exact diagonalization of very short chains, although recently there has 
been a significant numerical progress using time-dependent DMRG (Density 
Matrix Renormalization Group) procedure [58^ . 

Nontrivial exact solutions for time-dependent correlations do exist for free 
fermionic spin chains [59l [60l [611 [29l |62l [63j . due to the non-interacting nature 
of the excitations. For such chains, not only bulk regimes were investigated but 
also boundary effects [Ml [65l [Ml [SH [68] ■ On one hand, the z — z correlations are 
easily calculable due to their local expression in terms of the Fermi operators. 
They are basically fermion density correlation function. On the other hand, 
the X — X correlations are much more involved since in the Fermi representation 
one has to evaluate string operators. Nevertheless, as for the static correlators, 
one may use Wick's theorem to reduce them to the evaluation of a Pfaffian, or 
determinant, whose size is linearly increasing with i + j. 

For the A'^-sites free boundary isotropic XK-chain in a transverse field. 




J N-l ^ N 




the basic contractions at inverse temperature (3 are given by [64j: 

{Aj(t)Ai) = jpp^^^sinqj svnql ( cos Egt — i sin e^t tanh ^ 
{Aj(t)Bi) = sinqj singZ ( isine^t — cos£qttanh ^ 



(2.22) 



and the symmetry relations 



{B,{t)Bi) 
{Bi{t)A,) 



(2.23) 



where the excitation energies £q — J cos q—h with q — mi/ (A^+1) , n = 1, . . . , A^. 
With the help of these expressions, one is able to evaluate the desired time- 
dependent correlations, at least numerically. 
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The z — z correlator, namely (cr|(t)i7f ), is given in the thermodynamic Umit 
at infinite temperature T = oo by |69[ 164] 

(a|(i)af) [J,-i{Jt) i-iyjMJtf (2.24) 

where J„ is the Bessel function of the first kind. One has to notice that this 
result is field independent. The bulk behaviour is obtained by putting Z — s- cx) 
and keeping I — j finite. One has 

{a]it)an = JUiJt)''^h-\ (2.25) 

leading to a power law decay in time. For large time, the boundary effects lead 
to 

(a|(i)af) ^ ^ [sin^ {-^^ - y)] + ^^'^^^ 

with r = j — I the distance between the two sites. Therefore, the decay in time 
changes from t^^ to t^^ near the boundaries |64j . 

In the low-temperature limit (T = 0), we have a closed expression for h > J 
which is time and site independent |64j : 

K(^K> = 1 (2.27) 

revealing an ordered ground state. This can be related to the fact that the 
ground state corresponds to a completely filled energy-band, since all the energy 
excitations are negative. As already stated before, the calculation of the z — z 
correlation function is an easy task and we will not go on other models. 

In order to calculate the x — x time dependent correlation functions for such 
free-fermionic chains, one has to evaluate 

{cT^{t)af) = {Ai{t)B,{t)...Aj_At)B,^i{t)Aj{t)A,B,...Ai_,Bi_^Ai) . 

(2.28) 

In the high temperature limit (T — oo), the bulk correlation function of the 
isotropic XY chain in a transverse field h is given by the Gaussian behaviour 

[ZQlEileo]: 

4 

where 5ji is the Kronecker symbol. The boundary effects are hard to be taken 
into account due to the fact that the Pfaffian to be evaluated is a Toeplitz deter- 
minant [Tl] [59l [72] that can be treated only for large order [59] . Nevertheless for 
a vanishing transverse field, a conjecture was inferred from exact calculations 
for the boundary nearest sites i = l,2,...,5, claiming an asymptotic power-like 
decay [65| : 

{<J^it)af) - S,it-^/^-^'-^'^^^+^'^ . (2.30) 

At vanishing temperature and for h > J, the basic contractions can be 
evaluated in a closed form too. This leads to [73l [74l [64l [Ml 



(crj(i)crf ) = 5ji cos ht exp — (2.29) 



{<J^{t)af) = exp(-iM)exp[i(j - l)]J,-i{Jt) 



(2.31) 
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for the bulk correlation function. The Bessel function gives rise to a 
asymptotic behaviour. So moving from T = oo to T = there is a dramatic 
change in the time decay behaviour [64]. Finally, one is also able to take into 
account the free boundary effects. In this case, for large enough time, the 
behaviour changes from to t^^^^ ^64,- 

At finite temperature the asymptotic decay of the time- dependent correlators 
is exponential |75j : 

Kit)'Jt+n) « t'^'''-+''-^expf{n,t) (2.32) 

where f{n, t) is a negative monotonically decreasing function with increasing T. 
The preexponents i>± are known functions [75j of the field h, the temperature 
and the ratio n/t. When T oo, the function f{n,t) diverges logarithmically 
indicating the change of the decay shape from exponential to Gaussian |66| . 

To end this section, one may mention that the decay laws for the x — x 
time-dependent correlation functions of the anisotropic XY chain are basically 
the same as for the isotropic case, namely power law at zero temperature and 
Gaussian at infinite temperature [72j . 



Chapter 3 

Aperiodic modulations 



3.1 Definition and relevance criterion 

Since the discovery of quasicrystals in the middle of the eighties [76] , extensive 
studies on Ising quantum chains with quasiperiodic or aperiodically modulated 
couplings have been done in the nineties [III [78l [79l [53l [39] , see also [80| for a 
review. The interest in the field of critical phenomena in such aperiodic quantum 
chains lied in the fact that they offered intermediate situations between pure 
and random cases. 

Let us first define the aperiodic modulation itself. Most of the aperiodic 
sequences considered were generated via infiation rules by substitutions on a 
finite alphabet, such that A S{A), B S{B), where S{A){S{B)) is 
a finite word replacing the letter A{B). Starting from an initial letter and 
generating the substitution ad eternam, one obtains an infinite sequence of 
letters to whom a coupling sequence of the chain can be associated by the rules 
A —f Ja, B ~* Jbt ■■■ If one considers a two letter sequence then the cumulated 
deviation A(L) from the average coupling J at a length scale L is characterized 
by a wandering exponent lo such that 



where 5 is the strength of the aperiodic modulation and where the wandering 
exponent uj is obtained from the substitution matrix M whose elements {M)ij 
give the number of letters aj contained in the word S{ai) [811182] . The wandering 
exponent is given by 



where Ai and A2 are the largest and next largest eigenvalues of the subsitution 
matrix. Near the critical point of the pure system, the aperiodicity introduces 
a shift of the critical coupling 5t ^ ^'^-i ^ ^-1^(^-1) -^^jiicii j^as to be compared 



L 




(3.1) 
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to the distance to criticality t This leads to Luck's criterion 

j-^t-^ $ = l + j/(c^-l) (3.3) 

and in the Ising universahty class since v ~ I, one obtains 

$ = . (3.4) 

When Lu < 0, the perturbation is irrelevant and the system is in the Onsager 
universality class. On the other hand for ui > 0, the fluctuations are unbounded 
and the perturbation is relevant. For w = 0, that is in the marginal case, one 
may expect continuously varying exponents. 



3.2 Strong anisotropy, weak universality 

It is not here the purpose to give an exhaustive view on what was done in 
the context of aperiodic modulation of the Ising quantum chain, but rather to 
exemplify some (simple) aspects of it which have not really been noticed before 
(as far as my knowledge goes). 

The Ising quantum chain in a transverse field is defined by the hamiltonian: 

L-l L 
1=1 i=l 

where the cr's are the Pauli spin operators and Ji^hi are inhomogeneous cou- 
plings. In marginal aperiodic systems, for which hi = h and A,; = Ai?-^' with 
Ai = Ji/hi and /i = 0, 1 generating the aperiodicity, an anisotropic scaling was 
found. For such systems, the smallest excitations scale at the bulk critical point 
as A ~ with the size L of the chain. In the marginal aperiodic case, it was 
shown numerically that the anisotropy exponent z is continuously varying with 
the control parameter R [83l [53] 

2(7?) =x™^(i?)+a;„,(i?-i) > 1 , (3.6) 

where Xm^ (R) is the magnetic exponent associated to rus — (erf )0 The observed 
symmetry in the exchange i? <-> 1/i? in p.6p was demonstrated in [79] for 
aperiodic systems generated by inflation rules, using a generalization of an exact 
renormalization group method introduced first in [38] and applied to several 
aperiodic systems in |39| . We show here that this equation comes from a relation, 
valid for any distribution of couplings leading to anisotropic scaling, that rely 
the first gap Ai to the surface magnetization. 

Using a Jordan- Wigner transformation the hamiltonian (|3.5p can be 
rewritten in a quadratic form in fermion operators. It is then diagonalized by a 

^One may notice that relation l|3.6|l holds for the homogeneous system, R = 1, with 
^nig =1/2 and 2 = 1. 
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canonical transformation and reads 

L 



(3.7) 



9=1 



where 77^ and rjq are the fermionic creation and anihilation operators. The one 
fermion excitations Aq satisfy the foUowing set of equations: 

A,*,(i) = -h.'Pqii) - J,^q{l + l)A,$,(l) = - J._l*g(l - 1) - /l.,*,(l) (3.8) 

with the free boundary condition Jq — Jl = 0. The vectors $ and 5* are related 
to the coefficients of the canonical transformation and enter into the expressions 
of physical quantities. For example the surface magnetization, = (af), is 
simply given by the first component of $1 associated to the smallest excitation 
of the chain, Ai . 

Let us consider now a distribution of the couplings which leads to anisotropic 
scaling with a dynamical exponent 2; > 1, as it is the case for bulk marginal 
aperiodic modulation of the couplings [531 IMl IZSl ISHl • Then the bottom 
spectrum of the critical hamiltonian scales as Aq ~ in a finite size system. 
According to Ref . [THl [Ml [SS] , the asymptotic size dependence of Ai (L) is given 
by the expressions 



Ai(i)-(-l)^^ 



L-l 



L-l i 

2=1 j = l 



(3.9) 



ML) 



( ,)L '^iiL) -pr , 1 



L-l 



L-l i 



(3.10) 



which are valid at the critical point and in the ordered phase. Noting in p.9p 
that 



L-l 



UK' 

1=1 

equation p.9p becomes 



L-l i 

i=i j=i 



L-l 

HA. 



1+EriA,^ 

i=i j=i 



(3.11) 



L-l i 



i + ERA; 



i=i j=i 



Now multiplying both sides of equation (|3.12p with p.lOp leads to 



L-l i 



-1/2 



L-l i 



-1/2 



Ai(i) 



i+EnA.^ 



i+EnA.^ 



(3.12) 



(3.13) 
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which is symmetric under the exchange A <-> 1/A. One recognizes in this expres- 

'1/2 



sion the surface magnetization [42] ms{L, {Xi}) = 1 + J2f=i ^ 
the quantum chain, so that one finaUy obtains 



of 



Ai(L) ~ m,(L, {AJ)m,(L, {A^i}) . (3.14) 

This relation connects a bulk quantity, Ai, with surface quantities, namely 
ms{L, {Ai}) and the surface magnetization of the dual chain, ms{L, {A-^}). 

Consider now a deterministic distribution of the chain couplings, {A^}, with 
hi = h and A.; — JR^^/h with fi following some sequence of and 1. The 
critical coupling Ac follows from the relation limi_foo 11^=1 {Jk/huf'^^ = 1 27j 
and gives Ac — R"^"" with poo the asymptotic density of modified couphngs Ai?. 
The modulation of the couplings introduces a perturbation which can be either 
relevant, marginal or irrelevant. For the Ising quantum chain the fluctuations 
around the average couphng A at a length scale L A(L) = J2k=ii^k — \)^L'^ 
governs the relevance of the perturbation [20]. According to Luck's criterion, 
if w < the fluctuations are bounded and the system is in the Ising universal- 
ity class. On the other hand for w > 0, the fluctuations are unbounded and 
one has to distinguish two different situations. The evaluation of the surface 
magnetization is related to the sum 

A-2Ji?-2"^ where rij is the number of 
modified couplings, Ai?, at size j. At the critical point, the sum can be rewritten 
asymptotically as X^j^i R^^^^ ■ Now assume that the coefficient B is positif 
(if not the roles of i? > 1 and i? < 1 are reversed in the following discussion) . 
For i? > 1 the previous sum is absolutely convergent for L — > oo and leads to a 
finite surface magnetization with exponentially small corrections in a finite size 
system. On the other hand for i? < 1, the sum is diverging exponentially and 
the surface magnetization is governed by the dominant term exp(— 2i31ni?i"), 
so that 

m^(i.,i?) ~ exp(-A(i?)L'^) , (3.15) 

with A{R) > 0. So from equation (j3.14p . in both case (i? > 1 or i? < 1) the 
first gap Ai will show an essential singularity leading to z = c»: 

Ai(i?,i) - exp (-i(i?)i") , (3.16) 

with i(i?) = A{R) for i? < 1 and i(i?) = A(l/i?) for i? > 1 since from (|3ll| 
Ai(i?) = Ai(l/i?). 

In the marginal case, corresponding to w = 0, that is a logarithmic divergence 
of the fluctuations, rij ~ pooJ + C'lnj where C is some constant, it can be shown 
that the surface magnetization scales at the critical point as 

TO,(L,{A,}c)^L-^'"=(^) (3.17) 

with an exponent (i?) varying continuously with the control parameter R. 
In fact the sum can be evaluated at the critical point using rij ~ PocJ + Clnj. 
One obtains J^f^, A-2Ji?-2„, ^ j^f^^ r-^c in j _ jL ^^^-2C in r _ ^i-2Ci„i?_ 
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This expression is only valid in the weak perturbation regime for i? ~ 1, that is 
in first order in Ini?. For a stronger regime one has to retain higher terms in the 
Uj expression. At this order, the surface magnetic exponent is a;,„^(i?) ~ 1/2 — 
Clni?. One can remark that for a sequence like the period-doubling one [831153]. 
XmSR) = Xm,{^/R) which implies C — Q (one can test this numerically) and 
then the former calculation gives Xm XR) ^ l/2 + 0{\n^R). Finally, from ((XTi)) 
and p.l7p one obtains the relation 

z{R) = x„,^{R) + x,nAR~^) ■ (3.18) 

The anisotropy exponent z is then given by one surface magnetic exponent x^^ 
which is a function of the perturbation strength. The symmetry R <-> l/R of z 
is due to the self-duality of the Ising quantum chain which implies for all bulk 
quantities the relation Q{{Xi}) = Q({A^^}). 

For a symmetric distribution of couplings with respect to the center of the 
chain, leading to XmA^) = ^^m, (1/-R) (see the period-doubling case [83l [5311791 
[38l[39]), one observes a weak universality for the surface critical behaviour. The 
bottom of the spectrum scales anisotropically as A ~ ^ ^ t^"^ where t 
measures the deviation from the critical point and v = 1 is the exponent of the 
longitudinal correlation length So that from (|3.14p 

m,(t) ~ (ef^ . (3.19) 

From anisotropic scaling, one obtains for the critical dimension of the surface 
energy density the scaling relation Xe,. = z + 2xm^ ■ Using the symmetry of 

Xm, one has 

es - ■ (3.20) 

We see that we recover the homogeneous surface exponents Xm, — 1/2, x^, — 2 
when the deviation from the critical point is mesured with ~ Ai. 

The same weak universality seems to hold for the bulk quantities. In fact, 
it was shown in Ref.[52j that the bulk energy density scales as e ~ for 
marginal aperiodic modulation of the couplings. So again the pure energy den- 
sity exponent = 1 is recovered. Here we have investigated the behaviour 
of the mean bulk magnetization = ^/ L^^m{i) for the period-doubling se- 
quence from finite size scaling analysis. The magnetization is evaluated at the 
bulk critical point for sizes up to L = 1024. Numerically the profiles are well 
rescaled on the same curve with the exponent Xmb(.R) = z{R)/8 but we notice 
that as the size increases the profiles are more and more decorated with a grow- 
ing fiuctuation amplitude. This suggeres that the finite size behaviour of the 
mean critical magnetization is given by: 

mf,^L-^/8_ (L--)'^' . (3.21) 

On the basis of the numerical datas, the magnetization profile is compatible 
with the form 



m{l,L)^ |sin(^Z/L) ""-"i. [A + G{l/L)] , (3.22) 
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2.5 I ' ^s* 




0.0 ' ' < ' < ' < ' < ' 1 
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Figure 3.1: Rescaled magnetization profiles of the period-doubling chain with 
r = 5. The insert gives the corresponding magnetization profile. 

where ^ is a constant and G{x) is a kind of fractal Weierstrass function with 
zero mean value whose fourier momentum are given by the period-doubling 
cascade. The sinus term is very general for the profiles of the Ising quantum 
chains and is related to the geometry of the system. This can be demonstrated 
explicitely for the pure case [54j and was obtained numerically for random Ising 
systems [M]. The only difference here is that we have not only a pure constant 
in the front of it but in addition a fractal function of zero mean which controles 
the local fluctuations, due to the aperiodic distribution, around some average 
environnement. 

In conclusion the weak universality observed in this systems implies that the 
knowledge of the anisotropy exponent z is sufficient to determine the critical 
behaviour of the system. 



Chapter 4 

Disordered chains 



4.1 Random transverse Ising chain 
4.1.1 RG-equations 

Quenched disorder, i.e., time-independant disorder, has a strong influence on the 
nature of phase transitions and in particular on quantum phase transition. Many 
features of the randomness effects can be observed in one-dimensional systems, 
for which several exact results have been obtained and a general scenario (infinit 
randomness fixed point) has been developped [85]. In particular, since the work 
of Fisher [23], the Id random transverse Ising model has been the subject of 
much interest. The model is defined via the Haniiltonian: 

i=l 1=1 

where the exchange couplings, Ji and transverse fields, hi, are random vari- 
ables. In the case of homogeneous independently distributed couplings and 
fields. Fisher [53] was able to obtained many new results about static properties 
of the model. The method used is based on a real space renormalization group 
(decimation like) method [M], which becomes exact at large scale, i.e., suffi- 
ciently close to the critical point (for a general extensive review of the method, 
one can refer to [85]). Later on, these results have been checked numerically 
and results have been obtained on the dynamical properties at the critical point 
and outside, in the so called Griffiths phase [SllMlll71IH71[HilMl[nni[ni]- 

Let us sketch the main properties of the random transverse Ising model: 
For the inhomogeneous model, the phase transition is governed by the quantum 
control parameter 

6=[lnhU-[\iiJU (4.2) 

where the [.]av stands for an average over the quench disorder distribution. 
The system is in a ferromagnetic phase in the region S < 0, that is when the 
couplings are in average stronger than the transverse fields. Otherwise, it is in 
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a paramagnetic phase. The random critical point corresponds to S = and is 
governed by new critical exponents, that is a new universality class. 

The renormalization procedure used by Fisher is a decimation-like method. 
It consists of the successive elimination of the strongest coupling (bond or field), 
which sets the energy scale of the renormalization (the energy cut-off), f2 = 
max{Ji,hi}. If the largest coupling is a bond, that is Cl = Ji, then typically 
the randomly distributed neighboring fields hi and /i^+i are much more smaller 
than Ji, i.e., Jj ^ hi,hi+i. Thus, one can treat perturbatively the two spins 
connected via the coupling Jj and replace them by a single spin, with double 
momentum, in an effective field 



If the strongest coupling is a field, = //j. , then the connecting bonds Jk-i 
and Jk to the neighboring spins are typically much smaller, i.e., hk S> Jfe-i, Jk- 
Prom the magnetic point of view, this central spin is inactive, since due to the 
very strong transverse field, its susceptibility is very small and basically it gives 
no response to the total magnetization. Therefore, we can decimate out this 
spin and connect the remaining neighboring spins by an effective coupling 



The renormalization consists of following the evolution of the relevant distribu- 
tions during the lowering of the energy scale Q.. The locality of the procedure 
and the fact that the couplings are random independently distributed, leads to 
integro-differential equations for the physical distributions that can be solved 
exactly. 

Let us examplify the discussion by considering the distributions of fields and 
bonds, respectively P{h, Vl) and R{J^ il). Under the lowering of the energy scale 
from O to J7 — dVl, a fraction dil (P(fi, f2) + R{^1, Jl)) of spins are eliminated. 
The total number of bonds with value J at the new scale f2 — rffi is given by an 
equation of the type 

N{n - dn)R{j, n-dn) = N{n)R{j, o) -i- 7V(r2)Fi„(j, n) - N{n)F„^tiJ, ^) 

where, iV(f2) is the total number of spins at scale Q, and Fin^out{J,^) arc the 
flux of incoming and outcoming bonds with value J at the initial scale il during 
the decimation. With N{n - dVL)/N{n) = 1 - dn{P{n,n) + i?(0,fi)), and 
accounting for the decimation process, one arrives at 





hk 



R{j, n - dn) [1 - dn{p{n, n) + R{n, n))] 




(4.3) 
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From duality, it is easy to obtain a similar RG-equation for the field distribution 
just by interchanging fields and bonds and P and R. Expanding, R{J, ft — dfl), 
one arrives at the integro-differential equations 

^ = i?( J, n) {Pifi, n) ~ R{n, n)) 
-p{n,n) f dj' R{j',n)R{^n,n)^ , (4.4) 

and a similar equation for the field distribution 
dP 

— = p{h, n) {R{n, n) - p{n, n)) 
-R{n,n) [ dh' p{h',n)p{^n,i})^ . (4.5) 

Jh h n 

These are the basic equations, supplemented with the initial distributions, that 
have been solved by Fisher [53]. The results obtained are asymptotically exact 
as the line of fixed points is approched as ^ 0. We will not here continue into 
the details of the solution but rather we invite interested people to see ref. j85] 
for a recent extensive review on strong disorder RG methods. To mention only 
few of them, the magnetization critical exponent /3 is found to be /3 = f — r 
where r = ^^'^ is the golden mean and the dynamical exponent z = 1/(2|A|) 
where A is the asymmetry parameter related to the relative strengths of the 
couplings and transverse fields. At the critical point, A = and this leads to 
the log energy-scale | InSlj ~ for excitation energies on a finite system with 
size L. The average spin-spin correlation function G = [(cf crf_^,^)] involves the 
average correlation length which diverges at the critical point as ^ ~ S^'^ with 
V = 2 which differs from the typical value which is vtyp = 1- These differences 
are due to the non self-averageness of the problem. 



4.1.2 Surface behaviour 

The critical behaviour at the surface of the disordered chain can be obtained in 
a simple way, making a mapping to a random walk problem [171 ISO]- Indeed, 
fixing the end-spin = +, the surface magnetization is given by the formula 

m 

... -1/2 

hi 



I 



i=l j = l ^ - 



(4.6) 



which is of Kesten random variable type. In the thermodynaniical limit L — > oo, 
from the analysis of the surface magnetization, one can infer that the system is 
in an ordered phase with finite average surface magnetization for 



6 = [In /i]a« - [In J]a,, < 



(4.7) 
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where [.]av stands for an average over the disorder distribution. The system is 
in a paramagnetic phase for S > 0. 

To be more concrete, let's put ah the fields constant: hi — ho Vi and let 
the couplings be defined as 

J, = A'' (4.8) 

with A > and where the exponents 9i are symmetric independent random 
variables. Due to the symmetry of the random variables 9, the quantum control 
parameter is simply given by 

5 = \nho (4.9) 

so that the distance to the transition is controled only by the field value. The 
parameter A controls the strength of the disorder. For example, the binary 
distribution 

nie)^l6{e + i) + ^6(6-1) 

corresponds to couplings either J = A or J = 1/A. In the next subsection we 
will discuss the case where the distribution Tr{0) is very broad with power law 
tails corresponding to so called Levy flights or Riemann walk. But at the present 
let us consider only normal distributions (satisfying central limit theorems). 
The average surface magnetization, [ms]av(i), can be written as|42j 

(\ -1/2 
Vexp -21nA{5]")-5^i} (4.10) 
j=0 ) 

where the index n referes the different disorder realizations and with = (5/ In A 
and where the random sequence = 0^"^ + (f^^ + • • • + ^j"'' and 5q"^ = 
by convention. Since the random variables 0^ are independent identically dis- 
tributed variables we can interpret Sj as the displacement of a one dimensional 
random walker starting at the origin 5o = and arriving at a distance Sj af- 
ter J steps. At each step the walker perform a jump of length Q distributed 
according to the symmetric distribution. In the strong disorder limit, A 3> 1, 
the leading contribution to the average magnetization comes only from those 
walks whose mean velocity Sjj j is larger than the drift velocity 8^ up to time 
L, since then 1 + X^j^i exp(— 21nA {Sj — S^jY) = 0(1) and the corresponding 
contribution to the average magnetization is of order one. On the other hand 
when the walker mean velocity becomes smaller than , it gives an exponential 
contribution to the sum X]j'=i exp(— 2 InA {S'j — (5toj}) which is then dominated 
by the maximum negative displacement Sjmax — S^jmax^ so that the contri- 
bution to [m^Jav is exponentially small, of order O (e~'"^ {\Sj„^a:c-s^^jmax}\y 
means that the average surface magnetization [ms]av(-^) is proportional to the 
surviving probability after L steps of a walker which is absorbed if it crosses the 
6wj line. Stated differently, it means that the typical walks will give an expo- 
nentially small contribution to the average magnetization which is completely 
governed by rare events, that is non crossing walks. The problem of the surface 



4.1. RANDOM TRANSVERSE ISING CHAIN 



37 



magnetization is then mapped onto a random walk problem with an absorbing 
boundary. 

In order to obtain the survival probability, we use Sparre Andersen formula 

MM 

F{z,S)^exp{P{z,S)) (4.11) 

where F{z, 6) ~ X]n>o 5)z'^ is the generating function of the survival prob- 
abilities F{n,6) of the walker after n steps such that the displacement of the 
walker Sj > jS for all j < n and P{z,5) = ^^yii^i^^^)/''^)^^ where P{n,S) 
is the probability that at step n, Sn > nd. At the critical point, (5 = 0, since 
the distribution is symmetric we have P{n,0) ~ 1/2 for n ^ 1 so that from 
Eq. (|4.1ip we have F(n, 0) which implies 

K]av(L) - L-'/^ , (4.12) 

that is — 1/2. 

One is also interested in the typical behaviour at the critical point of the 
surface magnetization, m*^P(L) = exp([lnTOs(L)]av)- In the strong disorder 
regime, the leading contribution to the typical magnetization comes from those 
walks who visit the negative axis, that is, walks not contributing to the average 
magnetization. For such walks, Inm^ ~ — InA IS'jmaa;!, so that 

m*^P(L) = exp (- InA [|^,>„,.|]av) (4.13) 

with [|5'jma3;|]av the avcragc maximum displacement of the walker on the neg- 
ative axis after L steps. This average is actually given by the absolute mean 
displacement which scales at large L as L^^^ and which measures the fluctua- 
tions of the walk. So that the typical magnetization is 

m*^P(L) ~ exp(-const.Li/2) , (4.14) 

which means that the appropriate scaling variable is (In ms)/L^/^ and one can 
show that the distribution of nis is logarithmically broad. We close here the 
discussion of the surface critical properties of the (normal) random Ising quan- 
tum chain. One can see ref. [94j for exact results on the distribution functions 
of the surface magnetization. 

4.1.3 Levy flights 

In the present subsection, we study the influence of very broad disorder distri- 
butions on the critical properties of the random Ising quantum chain [5lJ . For 
that, we use a distribution tt{9) such that for large arguments it decreases with 
a power law tails 

7r{9) - \9\-^-" , |0| > 1 (4.15) 

where a is the so called Levy index. We will consider here the region a > 1 where 
the k-th moment of the distribution exists for k < a. More precisely, we consider 
the distribution Tr{9) = pa{l + 9)-'^-'' for 61 > and Tr{9) = qa{l + \9\)-^-°' 
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ioi 9 < 0, p + q = 1. In ref.[5T] we have also considered the discretized version 
(Riemann walk) of the above distribution. Using the continuous distribution 
defined above, the quantum control parameter is given hy S ~ In ho — {p — 
q) In A/ {a — 1). Taking ho = 1, leads to 5 — —{p — q) In A/(a — 1), which means 
that the asymmetry p — q drives the system outside the critical point <5 = 0. 

In the strong disorder limit A ^ oo, we have shown in the previous section 
that the average surface magnetization is given by the surviving probability of a 
random walk distributed according to tt{9) with an absorbing boundary moving 
with a constant velocity 6/\nA. So considering the following sum 

n 

Sn = Y.O,, (4.16) 

J=l 

in the large n limit, it exists a limit distribution p(u)du, in term of the variable, 
u = Sn/ln — Cn with the normalization 

^n^n^/" (4.17) 

giving the transverse fluctuation of the walk, if we interpret n = t as the (dis- 
crete) time and Sn=t as the position of the walker in the transverse direction. 
The second normalization is given by 

Cn = -n^-^'"5^ , (4.18) 

where with 5^ = ^(^) = ^/InA we define the bias of the walk. For a small 
asymmetry 5^ one gets from the combination in Eq. (|4.18p the scaling relation 
between time and bias as: 

t^\5^\-^^^\ ^[a)^^. (4.19) 
a — 1 

We use here the notation v for the exponent since remember that the final time 
t plays the role of the chain size L. For a symmetric distribution, when 7r(6') is 
an even function, thus p = q and = 0, we have for the limit distribution: 

p{u) = L„.o(u) = — / e^'="-l'=l"dfc , (4.20) 

which has an expansion around li = 



W.)^^£(-i)'|:^r(^) 



(4.21) 



fc=0 

and for large u it is asymptotically given by 

La,o{u) = -u-(i+")r(l + a) sin(7ra/2) , (4.22) 



where T{x) denotes the gamma function. 
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Consider next the surviving probability, Psurvit, 5^,), which is given by the 
fraction of those walks, which have not crossed the starting position until t = n, 
thus Si > for i = l,2,...,n. For a biased walk, with < \6w\ ^ 1, the 
asymptotic behaviour of Psurvin, 6^) is equivalent to that of a symmetric walk 
{Sw — 0) but with a moving adsorbing boundary site, which has a constant 
velocity of v — 6,^]. For this event, with Si > vi for i = 1, 2, . . . , n, the surviving 
probability is denoted by F(n,v), whereas the probability for Sn > vn, irre- 
spectively from the previous steps, is denoted by P{n, v) and the latter is given 
by: 

P{n,v) = / p{S,n)dS . (4.23) 



Between the generating functions: 

F{z,v) = ^^^(n,^;)z'^ 



PM . 5^^." (4.24) 



Pin, v) 

n>l 

there is a useful relation due to Sparre Andersen [HH [23] : 

Fiz,v) = e^p[Piz,v)] , (4.25) 

which has been used recently in Ref . [5S] . 

In the zero velocity case, v — 0, which is equivalent to the symmetric walk 
with Syj — 0, we have P(n, 0)=l/2. Consequently P(z, 0) = — iln(l — z) and 
F{z,0) = (1 — z)~^/^, from which one obtains for the final asymptotic result: 

Psurvit,0)= F{n,0)\^^^^t-\ ,9 = 1/2. (4.26) 

Note that the persistence exponent, -d = 1/2, is independent of the form of a 
symmetric probability distribution, tt{0), thus it does not depend on the Levy 
index, a. Thus, from the correspondence between the random walk problem 
and the random quantum chain, we obtain 

K(0,L)],, . (4.27) 

Thus the anomalous dimension xf^ = 1/2 of the average surface magnetization 
does not depend on the Levy index a and its value is the same as for the normal 
distribution (a > 2). 

For V = 6w > (paramagnetic case), i.e. when the allowed region of the 
particle shrinks in time the correction to P(n, 0) = 1/2 has the functional form, 
P{n,v) = l/2-g(c), with c = un^^i/". Evaluating Ea. (|4?23| with Eq. (|4?2T|l one 
gets in leading order, P{n, v) = ^ — cA{a) + O(c^), with A{a) = r(l + l/a)/7r. 
Then P{z,v) — P(z,0) ~ A(a)'y z"n^^/" is singular around z ^ 1^ as 

~ (1 — z)^'--^^^/°'\ consequently 

F{z, v)^{l- z)-i/2 cxp \~A{a)v{l - z)-^!-!/")] , (4.28) 
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in leading order and close to z = 1~ . Here the second factor gives the more 
singular contribution to the surviving probability, which is in an exponential 
form: 



r^^^exp -constJ^i^-i/" 



(4.29) 



Thus the average surface magnetization [ms{S,L)]^^ ~ Psurv{t,6w)\f^]^ has 
an exponentially decreasing behaviour as a function of the scaling variable 
Consequently the characteristic length-scale in the problem, the av- 
erage correlation length, ^, and the quantum control parameter, 6, close to the 
critical point are related as ^ ~ J-^^C") with 

^(«) = ^ ■ (4.30) 

Note that i^(a) is divergent as a 1+, which is a consequence of the fact that 
the first moment of the Levy distribution is also divergent in that limit. In 
the other limiting case, a — > 2~, we recover the result v — 2 for the normal 
distribution. 

For V = 6w < (ferromagnetic case), i.e. when the allowed region of the 
particle increases in time we consider the large \v\ limit and write Eq. (|4.23p 
with Eq.g^l) as P{n,v) ~ 1 - B{a)£-" + 0{£-^°') with B{a) = r(l + 
a) sin(7rQ!/2)/7ra. Then, in the large \v\ limit 

Piz,v) = -ln(l-z)-B(a)|w|-"^z"n-" , 

n>l 

where the second term is convergent even at z = 1. As a consequence the 
surviving probability remains finite as n —^ oo and we have the result, F(n, v) ~ 
1 — const|u|~", for \v\ ^ 1. For a small velocity, < |w| ^ 1, we can estimate 
F{n^ v) by the following reasoning. After n = Uc steps the distance of the 
adsorbing site from the starting point, ys = vric, will exceed the size of transverse 
fluctuations of the walk in Eq ((iTf|) . kr n^", with \v\-''^°''> . Then the 

— 1/2 

walker which has survived until ric-steps with a probability ofric , will survive 
in the following steps with probability 0{1). Consequently 

\im Psurv{n,S^)^ lim Fin,v)l^s n^t ^ Kl"^"^^^ ■ (4.31) 

So the average surface magnetization is finite in the ferromagnetic phase and 
for small \S\ it behaves as: 

hm [m,{6,L)l^ ^ \Sf% f3s = • (4.32) 



^In the Cauchy distribution case, corresponding to the limit a = 1, one can calculate ex- 
actly the surviving probability [95], and one finds F(t,Sw) ~ arctan i5,„ ^ which leads 
close to the critical point to F{t, Sw) ~ t~^^^ exp{— const. Sw Int) which has to be compared 
to the general a case. 
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Thus the scaUng relation, Ps = xf^^^u, is satisfied. 

To end with this section, we will put some words on the typical behaviour 
and distribution function of the surface magnetization. At the critical point, 
the typical samples, which are represented by absorbed random walks, have a 
vanishing surface magnetization in the thermodynamic limit. For a large but 
finite system of size L, the surface magnetization is dominated by the largest 
negative axis excursion of the walker (in the absorbing region) , that is by the 
fluctuation of the walk which scales as L^^", so that 

m*^P(L) - exp(-const.Li/") (4.33) 

and the appropriate scaling variable is Inrris/ L^^" . In the normal distribution 
limit a = 2, we recover the L^/^ dependence. In ref. ^Slj we have also studied 
numerically the bulk magnetic critical behaviour. We have shown that the criti- 
cal behaviour is controlled by a line of fixed points, where the critical exponents 
vary continuously with the Levy index a up to the normal limit a = 2, where 
they take their Fisher's value. 

4.2 Inhomogeneous disorder 
4.2.1 Definition 

Up to now we have only considered homogeneous disorder, that is site indepen- 
dent distributions. But in many physical situations the disorder is not homo- 
geneous. The inhomogeneity can be generated by the presence of a boundary 
(a surface) or an internal defect which may locally induce a perturbation on 
the distribution of the couplings and/or fields. Let us consider surface induced 
inhomogeneities which are characterized by a power-law variation in the prob- 
ability distribution: 7r/(J) — 7r(J) ~ and/or pi{h) — p(h) ^ for I ^ 1, 
such that the local control parameter S(l) varies as [SO] 

6{l) = 5-Ar''. (4.34) 

7r(J) and p{h) are the asymptotic, far from the surface, distributions. This 
choice for the functional form of the inhomogeneous disorder is motivated by 
the analogy it has with the so-called extended surface defect problem, intro- 
duced and studied by Hilhorst, van Leeuwen and others, in the two-dimensional 
classical Ising model |961 H?! HH [221 1100| . This type of inhomogeneity has been 
later studied for other models and different geometries. For a review see Ref. [19] . 
Such perturbation is expected to alter only the surface critical behaviour. Due 
to the asymptotic decay, the bulk remain unperturbed by the presence of the 
inhomogeneity. So that here, the bulk critical behaviour is in the Fisher univer- 
sality class. On the contrary, depending on the decay exponent k, one would 
expect a modification of the surface critical behaviour. We will illustrate that 
by considering the surface magnetic behaviour since, thanks to the random walk 
mapping already discussed, we are able to obtain analytical results jSOj . 
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4.2.2 Surface magnetic behaviour 

For the inhomogeneous disorder, the local quantum control parameter has a 
smooth position dependence which, at the bulk critical point, is given by S{1) ~ 
—Al^^ according to Eq. (|4.34[) . The corresponding random walk has a locally 
varying bias with the same type of asymptotic dependence, Sw{l) = —Ay^l^^. 
Consequently the average motion of the walker is parabolic: 

t . 

yp{t)^-J2sUl)^-^t^--. (4.35) 

J_ Hi 

1=1 

Under the change of variable y{t) — > y{t) — yp{t), the surviving probability 
of the inhomogeneously drifted walker is also the surviving probability of an 
unbiased walker, however with a time-dependent absorbing boundary condition 
at y{t) < -yp{t). 

The surviving probability of a random walker with time-dependent absorbing 
boundaries has already been studied in the mathematical [lOlj and physical 
|102|I103] literature. In a continuum approximation, it follows from the solution 
of the diffusion equation with appropriate boundary conditions, 

^^P{y,t) = D-^Piy,t), P[-yp{t),t]=0. (4.36) 

Here P{y, t) is the probability density for the position of the walker at time t so 
that the surviving probability is given by 

/oo 
P{y,t)dy. (4.37) 
-Vp{t) 

The behaviour of the surviving probability depends on the value of the decay 
exponent k. For k > 1/2, the drift of the absorbing boundary in Eq. (|4.35p is 
slower then the diffusive motion of the walker, typically given by 

y,{t) ^ {Dty/^ , (4.38) 

thus the surviving probability behaves as in the static case. When k < 1/2, the 
drift of the absorbing boundary is faster than the diffusive motion of the walker 
and leads to a new behaviour for the surviving probability. For A^ > 0, since 
the distance to the moving boundary grows in time, the surviving probability 
approaches a finite limit. On the contrary, for A^ < 0, the boundary moves 
towards the walker and the surviving probability decreases with a fast, stretched- 
exponential dependence on t. Finally, in the borderline case k — 1/2 where the 
drift of the boundary and the diffusive motion have the same dependence on t, 
like in the static case the surviving probability decays as a power, Psmvi^w, t) ^ 
^-e(A^) ^ however with a continuously varying critical exponent. 
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4.2.3 Critical behaviour with relevant inhomogeneity 

Let us first consider the probability distribution of Inrris on finite samples 
with length L at the bulk critical point 6 = 0. According to the arguments 
given above for the typical magnetization, [InTOgJav is expected to scale as 
[In J|j^j^(J///i;)]av when the perturbation tends to reduce the surface order, i.e., 
when A < in Eq. (|4.34|1 . Thus one obtains: 

L 

[Inmsjav ^ ^([In J;]av - [ln/li]av 

1=1 

L 

= A^r^'^AL^-^. (4.39) 
1=1 

The typical magnetization, defined through [In rrisjav = ln[ms]typ, has a stretched- 
exponential behaviour, 

[wsltyp ~ exp(const AL^^'^) ,A<0 . (4.40) 

For the probability distribution of In nig , Eq. (|4.39p suggests the following scaling 
form: 

(In m„L) = j^P^ [In m.Li-'^] . (4.41) 

which was successfully tested numerically jSO^ • 

The asymptotic behaviour of the scaling function P„i (u) as u — > depends 
on the sign of A. For enhanced surface couplings, A > 0, there is a nonvanishing 
surface magnetization at the bulk critical point as L — *■ cxd, thus the powers of 
L in Eq. ()4.41|) must cancel and Hm^^o Pm{v) ^ v~^. For reduced surface cou- 
plings, A > 0, limtj^o Pm{v) = 0, indicating a vanishing surface magnetization. 

Next we calculate the average behaviour of the surface magnetization, which 
is determined by the rare events with = 0{1). In the random walk picture, 
starting with A ^ A^ > 0, for A^ <C 1 one defines the time scale t* for which 
the parabolic and diffusive lengths in Eqs. (|4.35p and (|4.38p are of the same 

order, yp{t*) ~ yd{t*), such that t* ^ A^,^^^^ The surviving probability 

can be estimated by noticing that, if the walker is not absorbed up to t* , it 
will later survive with a finite probability. Thus PsurviA^ > 0) ^ (i*)^^/^ ^ 
Ai/^^-^"^ and the average surface magnetization has the same behaviour at the 
bulk critical point: 

[m,]av ~ ^1/(1-2''), K < i, < A < 1. (4.42) 

In the case of a shrinking interval between the walker and the absorbing wall, 
A^ < 0, the leading behaviour of the surviving probability can be estimated by 
looking for the fraction of walks with y(t) > —yp(t) which is given by Psurv(^) 
exp[— const A^ t^~2'*]. Thus, for reduced surface couplings, the average surface 
magnetization has the following finite-size behaviour at the critical point: 

[m,]av ~ exp (-const L^-^") , k < i, A < . (4.43) 
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One may notice the different size dependence of the typical and average surface 
magnetizations at criticaUty in Eqs. (|4.40p and (|4.43[) . respectively. 

To obtain the (^-dependence of the average surface magnetization in the 
thermodynamic limit, we first determine the typical size of the surface region 
Is which is affected by the inhomogeneity for 6 < and ^ < 0. It is given 
by the condition that quantum fluctuation 6) and inhomogeneity {Al~'^) 
contributions to the energy are of the same order, from which the relation Is 
follows. Inserting L ^ Is in Eq. g^S]), we obtain 

K]av ~ exp[-const \S\-^^-^^^/^], k < \,A<0. (4.44) 

In ref. [50], we have also investigated the behaviour of the first gap of the 
quantum chain and of the surface autocorrelation function. From the first gap 
study, the relation between time and length scales is found to be 

InT-.^i"". (4.45) 

In order to obtain an estimate for the average surface autocorrelation function 
one may notice that the disorder being strictly correlated along the time axis, 
a sample with a finite surface magnetization rris — 0(1) has also a nonvan- 
ishing surface autocorrelation function Gs(t) ~ m^. Since the fraction of rare 
events are the same for the two quantities, the scaling behaviour of [Gs(L, r)]av 
can be deduced from the corresponding relations for the average surface mag- 
netization. For enhanced surface couplings, yl > 0, according to Eq. (|4.42|) . 
limT^oo[Gs(T)]av ^ ['Tis(j4)]av ^ ^i/(i-2k) criticality. For reduced surface 
couplings, A < 0, the finite-size critical behaviour follows from Eq. (|4.43l) with 
hm7-^oo[Gs(i, T)]av ^ \ms{L,Af\aN ^ exp[— const L^^^"]. Now, using the 
scaling relation ()4.45p . we obtain 

(4.46) 

in the thermodynamic limit. 

4.2.4 Critical behaviour with marginal inhomogeneity 

In the marginal case k = 1/2, the differential equation with absorbing bound- 
ary conditions in Eq. (|4.36p can be solved |103] in terms of parabolic cylinder 
functions of order D^{x) [104| . For zero bulk bias, 6nj — 0, the surviving 
probability has the asymptotic dependence Pgurvit) ^ a-nd the exponent d 
is such that 

i^2^(-2A„) = 0, (4.47) 

taking D — 1/2 for the diffusion constant in Eq. (|4.36p . In the limiting cases 
the solution takes the form 




(4.48) 



exp 



-const (In''')' 
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for ^ 1 while for large values of one obtains asymptotically: 



e = J-A^exp[~2Al]. 




(4.49) 



Then the random walk correspondence leads to the scaling dimension of the 
average surface magnetization, 



where one has to identify the relation between the quantum chain perturbation 
amplitude A and the corresponding random walk drift amplitude A^. For the 
extreme binary distribution, one has A = A^,, but for other distributions the 
relation could be different. One has to take a reference point from which the 
identification is made or one may try to give a convincing enough argument 
[50] . Nevertheless, the surface magnetization exponent thus depends on the 
distribution, however only through the value of the parameter A^ entering in 
Eq. (|4.50p . the functional form remaining the same. 

Next we discuss the properties of the probability distribution of In mj . Re- 
peating the argument used above for relevant perturbations, we arrive to the 
scaling form of Eq. (I4.4ip . however with k = 1/2. We note that the same scaling 
form remains valid for irrelevant perturbations, with k > 1/2, but the scaling 
function Pm{v) has different limiting behaviours when w — s- in the two cases. 
While for K > 1/2 it approaches a constant, in the marginal case it goes to 
zero as lim^j^^o Pm{v) i)~2a:„(A)+i_ this way one obtains the proper scaling 
behaviour for the average surface magnetization. 

For the distribution of the energy gap at the critical point, the scaling rela- 
tion in Eq. (|4.45p remains valid, with k — 1/2. 

The asymptotic behaviour of the average surface autocorrelation function 
can be determined at criticality by scaling considerations like in a previous work 
for the homogeneous case [86]. Here we make use of the fact that, as already 
explained for relevant perturbations, the average autocorrelation function has 
the same scaling properties as the average surface magnetization. Under a 
scaling transformation when lengths are rescaled by a factor 6 > 1, such that 
I' = l/b, the average surface autocorrelation behaves as 



where we used Eq. (|4.45p to relate the time and length scales. Taking now 
b = (Inr)^ we obtain: 



(4.50) 



[G,(lnT)]av = 6-<"[G,(lnT/6i/2)] 



(4.51) 



[Gs{r) 



(Inr) 



(lnT)-''=(^). 



(4.52) 



av 



The last expression follows from the exponent relation in Eq. (|4.50p where the 
average correlation length exponent is v = 2. This scaling behaviour has also 
been verified by numerical calculations [50] . 
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Chapter 5 

Non-equilibrium behaviour 



5.1 Equations of motion 

In this chapter, we turn to the out of equihbrium dynamical behaviour of the 
free fermionic chains. The presentation is of review type and most of the re- 
sults exposed here were obtain by others. Nevertheless, in order to be as clear as 
possible, we have presented the topic in a rational way, using the available knowl- 
odge at the moment. Of course, the original references are cited within the text. 
Non-equilibrium properties of quantum systems have attracted a lot of interest 
through decades |105j . Precursor studies on free-fermionic spin chains were per- 
formed by Niemeijer [Mj and Tjion [106j at the end of the sixties. Soon after, 
Barouch, McCoy and Dresden |1071 1108j solved exactly the Liouville equation 
for the Xy-chain and computed the time dependent transverse magnetization 
for several non-equilibrium situations. They have shown the occurrence of an 
algebraic relaxation instead of exponential as predicted by Terwiel and Mazur 
|109j using a weak coupling limit. More recently, a special focus was pointed 
on non-equilibrium quantum steady states, which are driven by some currents 
|110l lllll [Il2j . These studies are motivated by the fact that quantum systems 
have a natural dynamics, given by the quantum Liouville equation |113| 

where p is the density operator, Ht the Hamiltonian and where C is the quantum 
Liouville super-operator acting on the vector space of linear operators. The 
expectation value of an operator Q at time t is given by 

{Q)it)=Tr{pit)Q} 

where we have assumed that the density operator is normalised to one. The 
traditional way to study non-equilibrium properties of a quantum system is to 
couple it with a bath that can be described itself quantum mechanically, for 
example an assembly of harmonic oscillators [114j . The total Hamiltonian is 
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then splitted into three different pieces: 

Ht = Hs + Hj + Hb 

where Hs is the system Hamiltonian, Hb the Hamiltonian of the bath and Hj 
stands for the interaction between the bath and system dynamical variables. 
Now, average system's quantities are obtained with the help of the reduced 
density operator, that is, density operator traced over the bath dynamical de- 
grees of freedom: 

Ps{t) = Tn{pm- 

The expectation value of a dynamical variable Q of the system is given by 

{Q){t)^Trs{Qpsit)} . 

For a complete review of this approach one can refer to Ref. [115j . 

Another, more simple, route to out-of-equilibrium problems, is to investigate 
the relaxation of a non-equilibrium initial state, in which a closed system has 
been prepared at time t — 0. For a pure initial state, j^"), the initial density 
operator is just the projector 

p = I . 

The time evolution of the state \'^) is thus simply given by the Schrodinger 
equation and is formally given by 

= exp(-ifft)|*) 

where H is the Hamiltonian of the closed system. 

In the framework of free fermionic models, the easiest way to handle the 
problem is to solve the Heisenberg equations of motion [59l [72] : 

^^X = i[H,X]. (5.1) 

Since all the spin variables of the quantum chains can be expressed in terms of 
the Clifford generators {F}, the first step is to solve the equation of motion for 
these operators. In order to avoid unnecessary minus signs, let us redefine the 
generators in the form of Ref. 59J : 

= (-i)"-'ri) 

( = (-im) (5.2) 

where the chain sites run from i = —N to i = N with the thermodynamic limit 
N oo implicitly taken. With the anisotropic Hamiltonian 

^ = E + J^'^^'^f+l + hot] 
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the equations of motion are 

'TT^n = J^Yn^i + ,PYn+l — hYn 

at 

= - J"X„+i - + hXn . (5.3) 

The sohitions are hnear combinations of the initial time operators {Pn}'- 

m 

— ^ ^ rn^— m(^) m/'m(^) (^■^) 



with 



and 



1 

— / dq exp(— img) cos(e„i) 

J-7V 



fmit) ^ — I dq exp(-img) cos(eqi) (5.5) 



/i™(i) = r dg exp(~imq)5^^^^ [j"e'« - /i + J^'e'^^] (5.6) 
where the excitation energies are 

{{r + jy) cos q-hf + [r - jyf sin^ J ^''^ . (5.7) 

In the special cases of the critical Ising chain, corresponding to — h — 1 
and jy — 0, and the isotropic XX-chain with — Jy = 1/2, one obtains 
closed analytical expressions for the basic time-dependent operators T2] . These 
closed forms permit to analyse exactly the long time behaviour of transverse 
magnetization profiles, end-to-end correlations [116[ 1117] . two-time transverse 
correlation functions and so on. In order to illustrate the procedure we concen- 
trate on the critical Ising chain following Ref. |59]. The isotropic chain can be 
treated on the same footings. At = h = 1 and jy = 0, the equations of 
motion simplify into 

^-Pn — Pn-l — Pn+l , (5-8) 

where we have used the operator notations P in order to have more compact 
expressions. By iteration, one obtains for the Ith order term 

^)'^"-=E(-l)'( 'l)Pn-l+2k, (5.9) 

k=0 



k 



where the ( ^ ) are the binomial coefficients. One can now sum up the Taylor 



series 



Eyr (d^j E P^-^M'i^), (5.10) 

(=0 ■ ^ ^ fc=-oo 
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where J„ is the Bessel function of integer order n. Finally one has 

oo 

P„(<)= PkJn-k{2t), (5.11) 

fc— — oo 

which gives the time dependence of the basic operators. Putting Ft{j) — Jj{2t), 
the time evolution of the operators P„ are expressed as a discrete convolution 
product: 

Pnit) ^J2Ftin-k)Pk = (Ft*P)(n) , (5.12) 

k 

where the kernel for the critical Ising chain is simply the integer Bessel function 
Jk(2t). This formula is exact for an infinite chain. In the case of a finite size 
system, one has to take care for site indices close to the boundaries. We will 
come back later on this point. Alternatively, one can express in closed form the 
functions / and h introduced previously. One has for the critical Ising chain 

/,n(t) = ^ / dge-™« cos (2t sin |) = J2„(2i) 
hmit) = ^f dge-'("-i/2),gi^ (2*«in|) = J2m-i{2t) (5.13) 

leading together with (|5.4[) to the time-behaviour (|5.1ip . 

Having solved explicitly the equations of motion for the Clifford operators, 
we can now express the time dependence of physical quantities in terms of the 
initial-time P„ operators. For example, the transverse magnetization at site n 
is given by 

<(i) = iP2n+i{t)P2n{t)iY,Ft{2n+l - k)Ft{2n ^ k')PkPk' (5.14) 

kk' 

which can be rewritten 

<W = Y.[Pt{2in-p))-F,{2{n-p) + l)Ft{2{n~p)~l)]a; 
p 

+ iJ2*Ft{2n+l-k)Fti2n-k')PkPk' (5.15) 

kk' 

where ^* restricts the summation over all k and fc' except fc' = fc± 1. The case 
k — k' , with = 1, gives zero since 

^ Pt(2n -fl - fc)Pt(2n - fc) = ^ Jk{2t)Jk+i{2t) = Ji(0) = Q . 
fc k 

The equation (|5.15p will be our starting point for the study of the transverse 
profile time-evolution with specific initial state. 



5.1. EQUATIONS OF MOTION 



51 



In the same way, the operator at time t is given by 

2n 

a:{t)^i~irl[J2P,^J,^,^{2t) (5.16) 

j=Q kj 

where we have labelled the first site of the chain by ^ = and taken a bulk site 
n, that is infinitely far away from the boundaries. In the string Yik , only 
those terms with all kj different will give a non vanishing contribution, since if 
at least two identical labels appear, we have a factor 

Jj_fe(2<) Ji_fe(2t) ^ Sij 

k 

vanishing for i ^ j. Finally one arrives at 

<W = (-i)" E PkoPk,-Pk,^det{a,,) , (5.17) 

fco<fcl<---<fe2Ti 

with 

= J,-k, (2<) < J < 2n . (5.18) 

The developments given so far are of course valid for infinite homogeneous 
chains. In the case of inhomogeneous chains, as for example in the presence of 
disorder |117j , the time dependent coefficients of the linear development of the 
operators Pn{t) are no longer given in terms of Bessel functions. Moreover, for 
finite homogeneous chains, they are also terms proportional not only to J„_fe(2i) 
but also to J„+fc(2i), where the last terms play a significant role for the near- 
boundary behaviour. In order to take into account these facts, we present 
now quickly the general expressions |118j valid for finite chains of the type 
(|3.5p . By introducing the Clifford FJj operators we arrived a.t H = (l/4)r^Tr 
with the 2L-component Clifford operator given by T^^ = (Pj, Fj, r|^) and 

= (r^,r^). The diagonalizatlon of the Hamiltonian can then be performed 
by the introduction of the diagonal Clifford generators 7^ = (7^,7^) related to 
the lattice one by F^^ = J2q 4>q{i^)l\ and = i^qip)^"^ with and ■0 defined 
previously as the eigenvector components of the matrix T given in (jl.lOp . It 
leads to 

9 

The time dependence of the diagonal operators is then simply given by 7^ [t) = 

Ul{t)lqUq{t) with 

U,{t) = cxp (^7^,^) = COS ^ + ^l^l sin '-f . (5.20) 
Utilising the fact that {7^,7^/} — 26ijdqq', we obtain 



2 



(5.21) 
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where we have defined the pseudo-scalar product as 

{C\D) = ^{C\D} (5.22) 

with {., .} the anticommutator. The time-dependent lattice Clifford generators, 
r^(t), can then be re-expressed in terms of the initial time operators F with the 
help of the inverse transforms 7^ — (f)q{k)Tl and 7^ = J2k ^qi^)^k- Finally, 
one obtains 

m)^J2(^l\rUt))Tl (5.23) 

k,i 

with components 

i^lKit)) = 5]0,(fc)0,(n)cose,t 

{rl\rl{t)) = {rl\rl{-t)) = -Y,Mk)M^)sme,t 

i^lWlit)) = Y.^g{k)Mn)cose,t. (5.24) 

These general expressions are exact for all finite size free boundaries free fermionic 
quantum chains. 

Formally, since (F^|F;) ~ SijS^i, the set {F^} forms an orthonormal basis 
of a 2L-dimensional linear vector space £ with inner product defined by (.|.) = 
.}. Hence, every vector X £ 8 has a unique expansion X = ^ ■ j.(FyX)F^. 
The string expression XiX2...Xn, with Xj € f , is a direct product vector of the 
space £1 ® £2 ® ■■■ ® £n which decomposition is 

xiX2...x„= (n\i^i>-(n:i^n)n\-n:- (5-25) 

ii .ki , . . . ,in .kn 

With the help of this time-development and the specific solutions and ip, we 
are able to analyse any non-homogeneous finite chain. One has to solve first 
the eigenvalue equation TVq = eqVq and re-inject the (j) and if) into the basic 
components given in (|5.24p and then use the general expression (|5.25p . 

5.2 Time-dependent behaviour 

5.2.1 Transverse magnetization 

In this section we consider initial pure states of the form 

I*) = |...a(fc)(7(fc + l)...) 

where cr(fc) is the value of the z-component spin at site k. Such initial states are 
physically relevant since they are quite easily accessible by the application of 
a strong modulated magnetic field in the desired direction. On the theoretical 
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side, these states permit to obtain exact solutions jll91 1120j . The question we 
ask is how the magnetization will relax as time evolves |1201 1121i I122[ 11231 1116j 
and how topological defects, such as droplets or kinks will spread out [llOilllll 
112311118] ? In order to give an answer, we compute the transverse magnetization 
profile 

(*krWI*>= E {r[\\Tf{t)){Tll\r}it)){^\-iTl\Tl\^). (5.26) 

Since in the z-state j'l') = • • • (K) |o'(fc)) \(j(k +1)) (g) . . . the only non-vanishing 
contributions come from terms — iF^r^ = (t|, the expectation value of the trans- 
verse magnetization at time t in such a state is given by: 

i^Wfitm = E [{Tl\rf{t)){rl\rjit)) - {Tl\Tfit)){rl\rjit))] a{k) . (5.27) 

k 

Clearly, for a translation invariant Hamiltonian this equation can be rewritten 
as a discrete convolution product |118| : 



m{l, t)^^Gt{l- k)a{k) = (G* a) {I) , (5.28) 

k 

with m(Z, t) = {^\af{t)\'i>). As it is seen from (|5.15[) for the critical infinite Ising 
chain the kernel Gt is given in terms of Bessel functions: 

Gt{l) = F^{21) - Ft{2l + l)Ft{2l ~ 1) (5.29) 

with Ft(n) = Jn(2t) as already defined. Due to the different asymptotic proper- 
ties of the Bessel functions one has to distinguish between the cases n/t — v > 1 
and u < 1. For v > 1, corresponding in (|5.28p to a distance n = I — k between 
sites I and k larger than we are in the acausal region since the elementary 
excitations, travelling with velocity equal to one by appropriate normalisation 
of the Hamiltonian ()3.5p . have no time to propagate from the initial position 
/ up to site k. This is exactly what is seen from the asymptotic behaviour 
of the Bessel function J„(t) which is vanishing exponentially as exp(— A(w)n) 
with X{v) > for n > t. So the behaviour of the local magnetization will be 
completely governed by the local environnement, since we have a compact sup- 
port kernel, insuring the existence of the convolution product. Inside the causal 
region, v < 1, with the help of the asymptotic [124j for v ^ I 



\ cos p J \ ■Kv tan f3 
where S = i/(tan/3 — /?) — 7r/4, one obtains in the continmmr limit 

Gtivt) = ^giv) (5.31) 

with 

9iv) = l H^-^r^' H<1 . (5.32) 
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In the continuum limit, the local magnetization m(n,t) — mt{v) is then given 
by the convolution product 

mt{v) = {g'kat){v) (5.33) 

with the initial state function crtiv) = a(tv). 

The same analysis for the XX-chain [118] leads for the Green function g{v) 

to 

9iv)^l -A^~-r"' H<1 , (5.34) 
^ [ \v\>l 

that is an inverse square-root behaviour. From these expressions one is able 
to evaluate the long time behaviour of the relaxation process from any initial 
2-state. With a homogeneous initial state, rn(0) = 1, one obtains for the XX- 
chain m{t) = to(0) = 1 as it should be, since the dynamics is conservative. 
On the contrary, in the Ising case, one has from (|5.28p together with (|5.32p 
m{t) = 1/2, that is in the long time regime half of the initial magnetization 
remains in the z-direction. A more careful analysiiS leads to [HSl [TT6] 

m{t) = i + 1 Ji(4t) (5.35) 

so that the final constant is in fact reached with a power law behaviour t^'^l'^ . 
One remarkable property of the relaxation of the magnetization of the Ising 
chain is that the remaining half magnetization has a conservative dynamics. 
For example, if we start with a droplet of L down spins inside an environment 
of up spins, in the long time limit t L, the magnetization profile is given by 



mtiv) = g{v) 



-5iv) 



1 2L , r 

--— Vl^. (5.36) 

Z TTt 



The excess magnetization m'^ = m — 1/2 is a scaling function and spreads out 
into the bulk of the up-spins without loosing any weight. In fact, one can write 
down a continuity equation dtm'^{n, t)-\-j{n, t)—j{n-~l, t) = where the current 
density is given in the continuum limit by [118] 

j(x,t) = |m=(x,t). (5.37) 

For a general initial z-state the Fourier transform of equation (|5.28p is 

mtiq) ^ at{q)g{q) (5.38) 
where the kernel in Fourier space is 

gig) = Jo(27rg) (5.39) 

for the XX chain and 

m = ^ (5.40) 



^The continuum analysis performed here is valid up to order t ^. 
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for the Ising one. By inverse Fourier transform one obtains the desired magne- 
tization profile. 

It is to be noted here that in the XX case an additional step-like structure 
arises on top of the scaling function [125j . It was found that this step- like 
structure is related to the flow of localized flipped spins. Indeed, each step width 
broadens as t^^"^ while the height decreases as t^-^^^, keeping the transported 
magnetization within each step equal to one unit. One may notice that the 
relaxation of spatially inhomogeneous initial states has been treated for several 
variants of the XY quantum chains, modulated, dimerised and so on |126j . 
A slowing down of the relaxation may occur for fermionic models that have 
a gaped excitation spectrum at some special points of the modulation wave 
vector |1211ll22[|126j . Such study was also extended, by time-dependent DMRG 
method, to the XXZ-chains [127] . Starting with a kink in the z-direction, it 
was shown that the magnetization transport remains balistic with a scaling 
function ^{n/t) for small 5^ 5^ interaction while it vanishes for strong enough 
interaction. Disordered Quantum chains have been investigated numerically 
in Ref. [117j . For the transverse magnetization, they show similar behaviour 
to the homogeneous chain, that is an algebraic decay in time. During the 
time evolution, the spatial correlations are building up and at long time they 
reach a size-dependent constant depending on the distance from criticality |117j . 
This behaviour is related to the distribution of rare samples that are strongly 
correlated due to large domains of strong couplings. 

5.2.2 Boundary effects 

For an open chain, one has to take care about boundary effects that could for a 
near-boundary spin modify the relaxation behaviour compared to the bulk one. 
We will illustrate this with the behaviour of the transverse magnetization in the 
case of the XX chain. Consider as an initial state a droplet of i-down spins 
(in the z direction) at the boundary of a semi-inflnite chain, the rest pointing 
in the opposite direction. The transverse magnetization at site I and at time t 
can be written 

L 

m{l,t) = l-2^Ft{l,k) (5.41) 

k=l 

with the Green function Ft(l,k), using the symmetry properties of the basic 
contractions, given by 

F,ii,k) = m\Tim' + \{Ti\Tfm^ 

= [Ji-kit)-i-lfJi+kit)r (5.42) 

which is exactly the form given in (|2.24p for the z — z correlator {af{t)a^) 
at infinite temperature. The reason for this coincidence lies in the fact that 
the expectations of the string operators P2i for the correlator at infinite 
temperature appear in the very same form in the z-state expectations of the 
transverse magnetization. So that, this is also true for the Ising chain. From 
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these coincidence, we see that the long time behaviour of the Green function 
is Ft ~ near the boundary. In fact, from the asymptotic analysis of (|2.24p 
we can calculate explicitly the long time behaviour of the magnetization. The 
droplet magnetization, M^{t) = (1/2) X]fc=i '^ki^) spreads in the bulk, for t ^ L, 
as 



M^(t) = 




(5.43) 



which as to be compared to |123j 

, , L 1 

M'{t) = ---— (5.44) 

Z TT r 

for a bulk droplet. 



5.2.3 Two-time functions 

Two-time functions (Q(ti)Q(i2)) are of primary importance in characterising 
non-equilibrium dynamics. In particular, they show the phenomenon of aging, 
that is, the dependence of the correlation functions on both times ti and t2 > ti, 
where ti is usually called the waiting time and specify the age of the system. 
This is in contrast to the equilibrium situation where the dependence is only 
on the time difference t2 — ti. Usually, at large waiting times, two distinct 
regimes develop: (i) ai short time differences the correlations are time transla- 
tion invariant whereas (ii) at long time differences, the relaxation is very slow 
and depends on the waiting time too. Aging was first considered in ultra-slow 
glassy dynamics [128) but it has been also investigated in simpler systems, ei- 
ther classical [129] as well as quantum [1301 11231 11311 1116j . The first attempt 
in this direction on homogeneous short range quantum spin chain was done in 
Ref. [123]. 

The two-time non-equilibrium correlation function we consider is defined by 



a,{h,t2) = {^\<7'^{h)a]{t2m , (5.45) 

with a z-initial state j^*) as already defined and t2 > ti. In terms of the time- 
independent basic operators, the two-spins product af {ti)aj {t2) takes the form 



aKti)a]it2)^- E c-;--;X(ti,t.)r-r-r-r- (5.46) 

ilki.. .Z4fc4 



where the coefficients of that development are given by applying formula (|5.25p 
with erf (ii)cr|(i2) = -T^{ti)T}{ti)T]{t2)T]{t2). Using Wick's theorem, the ex- 
pectation in the |^) state is given by 

6 

Q,{h,t2)^Y.^l\h,t2) + {^\at{t,)\^){^\a]{t2m , (5.47) 

fc=0 
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with a time-translation invariant element 

T^^tuh) =< A,Aj >'^^'^<B,Bj >*^-*i - < B,Aj >*2-*i< AiBj >*^^*i 

(5.48) 

and non-invariant terms 





= - < A,Aj >*^-*i [BBf^f^ 


(5.49) 




= < AB, >*^-*i [BA\l^f' 


(5.50) 




= < B,Bj >*^-*i [AA]^f' 


(5.51) 




= -< B,Aj [AB]l]f' 


(5.52) 




= -[AA^f^[BB]^f^ 


(5.53) 




= [AB]^f^[BA\^f\ 


(5.54) 



with 

\CD]^t = E^4< ^"^^ ^'^-^J - < "^"^^ ^"^1 
fc 

(5.55) 

where we have used the short notations of Ref. [1161 that is to recall, Ai — 
r], Bi = -irf and for the basic time contractions < XY >*= {X\Y{t)) = 
i{Xt,y(t)} that can be read from Eq. ([Oi)) . The elements Tl?{ti,t2) satisfy 
the set of symmetry relations 

Tl\t,M)^n\hM) 
TliUM) ^ -Ti\h,h) 

T;\t,,t2) ^ 'Ti\t2,h) ^ ■ > 

Tl\t,,t2)=n\t2,h) 

Ti^{t,M)^n\t2,h)- 

From the previous equations, together with the symmetry relations, one obtains 
for the connected symmetrised spin-spin correlation, 

a,(il,t2) = i(*|K(ti),a|(t2)}|*> - (*|<7f(ti)l*>(*k|(t2)|*) , (5.57) 

the simpler expression 

C^,{h,t2) = n'{tiM) - \AA^f-\BB^f- + \AB^f'\BA^f- , (5.58) 

which is valid for all free fermionic spin chains, prepared in an initial z-state. 
One may notice that for an initial x-state [116^ , the formula for the correlation 
function is the same, apart that one has to replace the initial definition of the 
[CD] [TTC]. 
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5.2.4 Critical Ising chain 

We concentrate now on the critical Ising chain, since at this point the con- 
tractions are simpler and one is able to give exact closed analytical expressions 

m- 

In particular, the two-time autocorrelation function, that is, i = j, takes 
a very simple form for a completely ordered z-initial state. Together with < 
AAj >*= (-1)*+V2(,_j)(2i) and < A,B, >*= i(-l)^+^+i J2(,_j)+i(2i) for the 
infinite chain, one obtains 

C{h,t2) = J^{2{t2 - ti)) - i[/(i2 + ti) - g{t2 - h)f (5.59) 

with f{x) = J2{2x) + Jo (2a;) and g{x) — J2{2x) — Jq{2x). The two-time auto- 
correlation can be rewritten 



C{hM)^4{'^{t2-h)) 



(5.60) 



The dependence of the correlation function on both ti and t2 reflects the non- 
equilibrium behaviour of the system. However, in the long-time regime 3> 1 
with small difference times t = t2 — ti <^ ti one recovers time translation invari- 
ance, that is, the two-time function depends only on the difference r and decays 
algebraically as t^'^. This power law decay subsists with the same exponent for 
any value of the transverse field h |116j . 

By Kubo's formula [105] . the linear response function Rl^-{ti,t2) at site j to 
an infinitesimal transverse field at site i is given by the commutator 

i(vl/|[af(ii),a|(t2)]|*) (5.61) 
and, on the same lines as before, it is expressed at the critical field as 

Rt;^{tiM)^Q],{.hM) + Ql{tiM) (5-62) 



with 



Q]Ati,t2) = J2r-l(2T)J2.+l(2(ti +t2))^^^ 

tl +12 

~ J2r+i(2r)J2.-i(2(ti +t2))-^^^ (5.63) 

tl + 12 



and a time translation invariant element 

Ql{h,t2) = QI{t) = -^J2.„i(2r)J2.+i(2T) (5.64) 

with r = j — i the distance between the two sites i and j. Translation invariance 
in time is restored at small difference times, 1 ^ r = t2 — ti <C ii, since the 
dominant term is then Qf,, (t) and it leads to a t^'^ behaviour. 
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5.2.5 XX chain 

For the XX-chain, the dynamics is conservative, that is the total magnetization 
in the z direction is conserved during time evolution. Therefore, in order to 
consider a relaxation process one cannot use the same procedure as in the Ising 
case. Instead of an infinite chain fully magnetized in the up-direction, we take 
as the initial state of the infinite chain a droplet of L-spins pointing upward in z- 
direction, from site labels i = Otoi = i — 1, the rest of the spins pointing in the 
opposite direction |123j . This can be assimilated to a system of size L interacting 
through its boundaries with a bath, where it can dissipates magnetization. The 
final stationary state reached by the system is the completely magnetized state, 
with all the system's spins down, which is not the equilibrium state. 

The two-time connected autocorrelation function of the droplet is defined by 

C{h,t2) = {^\M%h)A'r-{t2m - {^\M'{t,)\^){^\M%t2m (5.65) 

with M^{t) = \ Y^^=o '^i^'^) the total transverse magnetization of the droplet 
subsystem. The 1/2 factor has been introduced in order to follow Ref. |123| . 
The Heisenberg equations are easily solved in this case and one obtains for the 
lattice fermions 

Ck{t) = e-'"* '^'^'^ Jk-r{t)Cr . (5.66) 
r 

Introducing these solutions inside (|5.65p . using Wick theorem and the initial 
droplet state, one obtains 

L-l L-l 



C{h,t2)= Jk-l{t2 - tl)Y Jk-r{t2)Jl-r{tl) 

k,l=0 r=0 

L-l 

- J2 Jk-r{t2)Jk-s{t2)Jl-r{tl)Jl-s{tl) . (5.67) 
fcj,r,s=0 

The long-time behaviour of that expression is obtained with the help of the 
asymptotic expansion of the Bessel function J„(t) ~ •\/2/7rfcos(t — n7r/2 — 7r/4). 
For ^1,^2 3> L, one obtains the asymptotic form 

CitiM)^ — . . (5.68) 

In the finite magnetization density regime 1123', that is for times smaller than 
L^, the second term dominates and one has finally 

(5(ii,t2) = -^^^tr'/(|) i«^l,^2«i^ (5.69) 

with f{x) ^ x^^ . One may notice that this scaling form is specific of aging 
at a critical point where critical coarsening takes place |132[ 11331 1134] . This 
means that aging occurs in the ordered isotropic XX-chmxi with deterministic 
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quantum dynamics [123j . In fact, to conclude positively about the occurrence of 
aging in such system, Schiitz and Trimper |123| have restricted their definition 
of aging to the phenomenon of increasing relaxation times, i. e. the longer one 
waits, the slower the relaxation becomes. For that purpose they considered the 
ratio 

R{ti,t2) = C{hM)/C{hM) = h/t2 , (5.70) 

that is a kind of normalized autocorrelation function with respect to the waiting 
time ti , which clearly shows that aging occurs in this homogeneous deterministic 
quantum system. 



Discussion and summary 



We have presented in this review the general "ahiiost" canonical approach to 
free fermionic quantum spin chains. Almost canonical in the sense that, several 
variants of notations and approaches appear all along the enormous literature 
devoted to this topic, and one is sometimes confused. However, we tried to give 
a systematic analysis and emphasise the basic technical tools, such as the Wick 
expansion, pair contractions and operator time development. 

We have presented very few results on the well known static properties of 
these spin chains but merely insisted on the way to calculate them, in this 
sense it is a pedagogical attempt. Moreover, we have shown explicitly how 
one can reconstruct the zero temperature quantum phase diagram from the 
study of the boundary magnetization. Since it is very easy to compute the first 
site magnetization [321 US] , one has paid a very low price for the knowledge 
of the bulk behaviour, compared to standard approaches [28]. Of course, the 
present analysis was a posteriori and, as everyone knows, it is an easy task and 
sometimes pathetic to recover old results by new approaches. However, one can 
imagine to use this bias to tackle really new systems, so to say to infer bulk 
behaviour from surface properties. 

Results on equilibrium time-dependent correlation functions are quickly re- 
viewed and I apologise for important works that have not been cited here for 
various reasons, the main being fortuitous omission or lack of knowledge. One 
should not worry too much about it, since the science is, or should be, a public 
good. Phrased differently, intellectual property is still a property. To summarise 
on the spin-spin time dependent correlations, the calculation of the transverse 
functions is an easy task, since they have a local expression in terms of the 
non- interacting fermions. They are basically fermion density correlation func- 
tions and show up power law decay in time. The x spin component correlations 
are much more difficult to evaluate, since in this case the Jordan- Wigner trans- 
formation leads to string Fermi operators. Nevertheless, using Wick's theorem 
their asymptotic behaviour can be extracted. At vanishing temperature, they 
decay with a power law whereas at infinite temperature they have a Gaussian 
shape. In the intermediate temperature regime, the decay is exponential with 
a power law prefactor. 

After the generally known properties of the homogeneous system, we have 
presented some results on aperiodically modulated systems and random ones. 
In the aperiodic case, we have briefly given a relevance criterion for such pertur- 
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bation and emphasized the anisotropic scahng behaviour found in such systems 
as well as the weak universality observed for some of them. In the random case, 
in a first part we summirized RG Fisher's results for the homogeneous normally 
distributed random Ising quantum chain. In the remaining part, we have pre- 
sented some analytical results obtained in cases where the disorder is no more 
homogeneous but presents a power law spatial variation or, in a second case, 
where the disorder is no more normally distributed but rather has long tails of 
Levy type. 

Finally, we dealt with some of the non-equilibrium properties of these sys- 
tems. We focused our attention, after solving the Heisenberg equations of mo- 
tions for the Fermi operators, on the relaxation of the transverse magnetization 
in the critical Ising chain and the isotropic XY chain. For some particular initial 
states, magnetized in a given direction, we derived analytically a general expres- 
sion that can be used to calculate at later times the transverse magnetization 
profile. In fact, the long time relaxation of the magnetization is expressed as a 
convolution product of the initial state with a response kernel obtained analyt- 
ically for both XX and Ising chain. It is seen in particular that the Ising chain 
exhibits some similarities with the conserved dynamics XX chain. That is, af- 
ter a transient regime, the relaxation process in the Ising case is conservative 
too. Two-time functions are also discussed in the light of aging phenomena. It 
is shown how, starting with some non-equilibrium initial state, the correlations 
depend explicitly on both the waiting time, characterising the age of the system, 
and the later measurement time. Although aging was first considered in com- 
plex systems such as structural glasses or spin glasses, it is in fact also present 
in simple homogeneous systems with a completely deterministic dynamics. This 
final result is in fact very recent and it has been demonstrated in classical (129| 
and quantum systems |123j . 
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